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ABSTRACT 


This  project  quantifies  several  aspects  of  a  new 
multiplexed  imaging  technique  proposed  by  D.S.  Davis.  The 
novel  approach  of  this  technique  involves  the  use  of  encoding 
masks  derived  from  a  basis  set  of  two-dimensional  Walsh 
functions.  There  were  two  distinct  problems  addressed  by  this 
thesis  research.  First,  a  study  of  computer  simulated 
diffraction  patterns  of  the  photon  flux  through  these  encoding 
masks  yielded  design  constraints  to  be  incorporated  into  a 
prototype  system.  These  constraints  were  expressed  in  a 
simple  mathematical  relation  in  terms  of  wavelength, 
diffraction  angle,  and  spatial  frequency.  A  second  problem 
addressed  the  minimum  spatial  resolution  required  for  pattern 
recognition.  The  conclusion  reached  was  that  the  minimiun 
number  of  resolution  elements  necessary  for  pattern 
recognition  is  64  in  each  direction.  This  determination  also 
fixed  the  minimum  size  of  the  basis  set  of  two-dimensional 
Walsh  functions  required  for  multiplexing,  as  well  as  the 
number  of  pixels  required  to  display  the  image. 
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Z .  IHTRODDCTIOII 


A.  MOTIVATION 

This  thesis  is  part  of  a  much  larger  research  and 
development  effort  being  conducted  by  the  author's  thesis 
advisor  [Ref.  l].  The  ultimate  goal  of  this  larger  project  is 
to  develop  and  evaluate  a  new  type  of  instrument  for  infrared 
imaging  and  imaging  spectroscopy.  The  key  feature  of  this  new 
technique  is  that  it  will  multiplex  images,  by  means  of  focal 
plane  encoding  masks,  onto  a  single  detector,  or  at  most  a 
pair  of  detectors.  The  viability  of  this  method  was 
demonstrated  in  the  proof -of -concept  thesis  research  conducted 
by  Capt  R.H.  McKenzie,  III,  USMC  [Ref.  2].  The  motives  for 
this  thesis  project  have  been  twofold.  First,  since  the 
masking  technique  will  be  employed  at  long  wavelengths,  a 
(quantitative  prediction  of  the  diffraction  behavior  is  needed 
before  mask  design  is  undertaken.  This  project  successfully 
addresses  that  problem.  The  second  motive  for  this  research 
project  is  a  direct  consequence  of  the  multiplexing  techni(que 
itself.  The  spatial  resolution  of  a  multiplexed  image  will  be 
inversely  proportional  to  the  number  of  encoding  masks 
incorporated  into  a  working  instriunent  [Ref.  1].  This 
project,  therefore,  successfully  provides  initial  estimates 
for  the  number  of  encoding  masks  needed  to  produce 
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recognizable  multiplexed  images  from  a  variety  of  commonplace 
objects  or  targets. 

B.  BACKORODMD 

Terrestrial  objects  emit  most  of  their  thermal  radiation 
at  infrared  (IR)  wavelengths.  Many  interesting  objects  are  in 
approximate  thermal  equilibrium  with  their  surroundings,  which 
have  temperatures  in  the  275  -  310  K  range.  Assuming  that 
these  sources  radiate  approximately  as  blackbodies  [Ref.  3], 
the  bulk  of  their  radiance  lies  within  the  5  -  25  nm 
wavelength  range.  Of  course,  warmer,  man-made  objects  will 
radiate  at  correspondingly  shorter  wavelengths.  It  is 
desirable  to  utilize  this  abundance  of  radiated  energy  to 
detect  and  to  identify  objects.  Infrared  images  provide 
infonaation  about  the  size,  shape,  and  temperature  of  the 
radiating  sources.  Infrared  spectra  reveal  the  chemical 
composition  and  thermodynamic  information. 

D.S.  Davis,  advisor  on  this  thesis  project,  has  been 
involved  with  state-of-the-art  spectroscopic  IR 
instrumentation  for  many  years.  Historically,  his  research 
interests  have  involved  the  development  and  use  of  unique, 
high  resolution,  fully  multiplexed,  infrared  Fourier  transform 
spectrometers.  His  interest  in  IR  multiplexing  techniques 
inspired  the  concept  of  multiplexed  imaging,  which  is  one  of 
the  goals  of  the  larger  research  project  of  which  this  thesis 
is  a  part.  [Ref.  1] 
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The  purpose  of  infrared  imaging  is  to  measure  the  image's 
irradiance  as  a  function  of  two  independent  spatial  variables, 
such  as  conventional  Cartesian  (x,y)  coordinates. 
Spectroscopic  imaging  introduces  a  third  independent  variable, 
such  as  wavelength  or  freguency.  Traditional  IR  imaging 
devices  and  imaging  spectrometers  typically  scan  over  one  or 
more  of  the  independent  variables  which  describe  image 
irradiance.  There  is,  however,  an  inherent  inefficiency 
associated  with  such  scanning  techniques.  This  inefficiency 
stems  from  the  fact  that,  by  observing  only  one  small  element 
of  the  image  field  (and/or  one  spectral  element)  at  a  time, 
the  instrument  ignores  most  of  the  energy  from  the  scene. 
This  can  have  severe  consequences  in  the  infrared,  because 
available  photons  from  a  target  are  precious,  and  they 
typically  comprise  only  a  small  fraction  of  the  entire, 
background-dominated,  IR  irradiance.  The  goal  of  the  proposed 
technique  is  to  develop  a  more  efficient  generation  of 
instruments  which  will  apply  the  proven  advantages  of  infrared 
Fourier  transform  spectroscopy  (FTS)  to  infrared  imaging  and 
infrared  spectroscopy.  [Ref.  1] 

C.  STATEMENT  OF  THESIS  PROBLEM 

This  thesis  research  addresses  two  fundamental  (questions 
which  must  be  answered  before  the  design  of  a  prototype 
multiplexed  imaging  instrument  is  complete.  First,  since  the 
instrument  will  operate  at  comparatively  long  wavelengths. 


3 


there  is  concern  that  diffraction  effects  produced  by  the 
encoding  masks  will  cause  substantial  beam  spreading  and 
degradation  of  image  resolution.  Therefore,  the  initial 
portion  of  this  thesis  project  quantifies  the  diffraction 
behavior  of  the  masks,  so  that  this  behavior  can  be  included 
in  subsequent  instrumental  design  considerations.  Second,  the 
image  spatial  resolution  required  to  produce  identifiable 
images  of  ordinary  terrestrial  objects  must  be  determined. 
When  the  proposed  encoding  scheme  operates  at  short 
wavelengths  (far  removed  from  the  diffraction  limits 
considered  in  the  first  part  of  the  research) ,  it  is 
conceptually  capable  of  extremely  high  spatial  resolution.  In 
fact,  the  resolution  will  be  inversely  proportional  to  the 
number  of  encoding  masks  used.  Determination  of  the  minimum 
useful  number  of  masks  is  important  so  that  needless  effort 
and  resources  are  not  wasted  on  an  over-designed  prototype 
instrument . 
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II.  OVERVIEW  OF  MULTIPLEXED  IMA6IM0 


This  section  provides  a  synopsis  of  the  multiplexed 
imaging  technique.  Although  the  ultimate  goal  of  the  entire 
project  is  to  develop  a  multiplexed  imaging  spectrometer,  the 
thesis  research  at  hand  considers  only  multiplexed  imaging. 
Therefore,  the  spectroscopic  aspects  of  the  larger  project  are 
ignored . 

A.  MULTIPLEXING  CONCEPTS 

The  key  feature  of  optical  multiplexing,  and  the  one  that 
produces  its  inherent  efficiency,  is  that  all  of  the  photons 
from  the  target  are  detected  all  of  the  time.  This  is 
achieved  in  a  multiplexed  scanning  instriunent  by  detecting 
linear  combinations  of  the  signals  of  interest,  rather  than  by 
scanning  through  sets  of  the  individual  signals  themselves. 
For  instance,  consider  the  operation  of  a  raster  scanning 
imager,  such  as  a  FLIP.  Let  0,  be  the  radiative  flux  through 
the  ith  pixel  of  the  image,  and  be  the  instantaneous  power 
sensed  by  the  detector .  Assume  that  there  are  no 
inefficiencies  or  losses  associated  with  the  detection 
process.  Then,  for  N  pixels  which  comprise  the  entire  image. 


P,-  -  0,.  ,  {i  -  1,...,N} 


(1) 
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While  the  device  observes  the  ith  pixel,  it  ignores  the  other 
N  -  1  pixels. 

A  multiplexed  sensor  detects  a  linear  combination  of  the 
pixel  fluxes 

N 

,  {i  -  , 

j=1 

(2) 

Where  the  are  weighting  coefficients.  Therefore,  one  scans 
through  N  /x's  rather  than  N  ^'s.  The  multiplexing  weights  fi.. 
are  not  chosen  arbitrarily;  there  are  several  constraints 
which  they  must  meet  [Ref.  l].  First,  equation  (2)  must  be 
invertible,  so  that  the  individual  pixel  fluxes  ip.  can  be 
recovered  from  the  measurements  of  the  P^.  Second,  the 
need  to  permit  as  much  radiation  as  possible  to  pass  to  the 
detector,  in  order  to  maximize  the  overall  optical  efficiency. 
The  inherent  inefficiency  of  traditional  raster  scan  devices 
is  due  to  the  fact  that  their 


Mfj  -  ^ij  f 


where  is  the  usual  Kronecker  delta  symbol.  Third,  the 
must  be  capable  of  being  implemented  in  a  practical  optical 
design  which  is  suitable  for  a  working  instrument. 

For  example,  consider  the  process  by  which  a  Fourier 


transform  spectrometer  (FTS)  encodes  spectral  intensities 
[Ref.  1].  In  this  instance,  the  <p.  represent  spectral  fluxes 
at  N  different  frequencies  or  wavelengths.  The  /x-j  are 
discrete  samples  of  the  familiar  orthogonal  sinusoidal 
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functions  of  Fourier  analysis.  The  are  the  discrete 
samples  of  the  recorded  interferogram.  Inversion  of  equation 
(2)  to  yield  the  spectral  fluxes  is  accomplished  by  applying 
a  discrete  Fourier  transform  to  the  P,. 

B.  SEQDEMCY  THEORY  RMD  MULTIPLEXED  IMAGING 

The  family  of  functions  chosen  to  implement  the  encoding 
scheme  for  this  multiplexed  imaging  project  are  the  Walsh 
functions  from  sequency  theory  [Refs.  1,4]  rather  than  the 
sinusoidal  functions  associated  with  FTS  and  Fourier  analysis. 
The  Walsh  functions  are  similar  to  the  sinusoids  in  that  they 
constitute  complete,  orthogonal  basis  sets.  In  other 
respects,  however,  they  are  radically  different.  The  most 
striking  differences  are  that  they  are  discontinuous,  and  that 
they  assume  only  one  of  two  possible  values,  ±  i.  Another 
difference  is  that  the  Walsh  functions  are  not  necessarily 
periodic.  Furthermore,  a  complete  basis  set  of  Walsh 
functions  must  contain  N  «  2”  members,  where  m  is  a 
nonnegative  integer.  Detailed  discussions  of  these  and  other 
properties  may  be  found  in  Beauchamp  [Ref.  4]. 

The  concept  of  a  Walsh  function's  sequency  is  fundamental 
to  understanding  this  thesis  research.  Its  definition  is 
straightforward.  The  domain  over  which  a  basis  set  of  Walsh 
functions  is  defined  is  partitioned  into  N  subdomains  of  equal 
length,  where  N  is  the  number  of  functions  in  the  complete 
basis  set.  The  Walsh  functions  then  assume  constant  values. 
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either  +1  or  -1,  across  each  subdomain,  according  to  the  rules 
which  dictate  their  generation  [Ref.  4].  At  the  boundaries 
between  adjacent  subdomains,  a  function  may  switch  values, 
thereby  crossing  zero.  The  total  number  of  such  zero 
crossings  that  occur  across  the  entire  domain  of  the 
function's  definition  is  that  function's  seguency.  For  a 
complete  Walsh  basis  set,  the  corresponding  member  function's 
seguencies  range  from  zero  to  N  -  1.  N  is  referred  to  as  the 
order  of  the  basis  set. 

C.  THE  NEW  IMAGE  ENCODING  TECHNIQUE 

The  new  multiplexed  image  encoding  technigue  will  utilize 
Walsh  functions  to  generate  optical  encoding  masks.  These 
masks  will,  in  turn,  be  used  to  modulate  the  flux  distribution 
in  an  image,  in  analogy  with  eguation  (2) .  An  illustration  of 
how  this  might  be  done  is  provided  by  the  following  example. 
Suppose  that  one  desires  to  determine  the  radiant  fluxes  0j 
through  different  contiguous  regions  of  an  image  field. 
Consider  the  complete  Walsh  basis  of  order  4.  Construct  a 
one-by-four  row  matrix  for  each  basis  function  in  the  set, 
such  that  each  function's  subdomain  structure  corresponds  to 
one  entry,  or  matrix  element: 
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Seauencv  Row  Matrix 


0 

(+1 

+1 

+1 

+1) 

1 

(-1 

-1 

+1 

+1) 

2 

(-1 

+1 

+1 

-1) 

3 

(+1 

-1 

+1 

-1) 

The  reader  may  easily  verify  that  the  rows  are  mutually 
orthogonal;  in  fact,  so  are  the  columns.  Hence,  the  rows 
could  be  combined  into  a  single  four-by-four  orthogonal  Walsh 
encoding  matrix,  such  that 


+1  +1  +1  +r 
-1  -1  +1  +1 
-1  +1  +1  -1 
,+i  -1  +1  -1, 


(4) 


The  practical  optical  implementation  of  this  matrix  would 
involve  construction  of  four  masks,  each  of  which  is  portioned 
into  four  regions  where  boundaries  match  those  of  the  image 
regions  to  be  multiplexed.  A  possible  example  of  such  a 
portioning  is  shown  in  Figure  1.  Each  mask  in  Figure  1 
corresponds  to  one  row  of  matrix  (4) .  Dark  regions  imply  a 
matrix  element  of  and  signify  that  the  mask  is  optically 
transmitting  across  those  pixels.  Light  regions  correspond  to 
matrix  elements  of  -l  and  denote  that  those  pixels  are 
optically  reflective.  Light  that  is  transmitted  is  collected 
and  sent  to  one  detector;  that  which  is  reflected  is  sent  to 
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a  second  detector.  The  second  detector's  output  is  subtracted 
from  that  of  the  first  detector.  As  the  encoding  instrument 
cycles  through  the  four  masks,  the  differenced  detector  signal 
will  be  proportional  to  the  Pj  in  equation  (2) .  Collecting 
the  four  P/s  into  a  single  column  vector  and  multiplying  it 
by  the  inverse  of  would  yield  the  desired  0j's. 


Figure  i  Two*>dimensional  Seguency  -Ordered  Walsh  Masks  for 

Four  Pixels 

The  actual  encoding  scheme  in  a  prototype  multiplexed 
imaging  system  will  not  make  explicit  use  of  the  two- 
dimensional  mask  configurations  like  those  in  Figure  1. 
Rather,  the  prototype  will  utilize  two  essentially  one¬ 
dimensional  mask  families,  and  it  will  exploit  a  very  powerful 
mathematical  feature  of  the  Walsh  functions,  the  Kronecker 
product  [Ref.  4].  Via  this  approach,  a  two-dimensional  mask 
basis  may  be  synthesized  from  two  lower-order  one  dimensional 
bases.  This  is  again  illustrated  by  example.  Consider  the 
Walsh  basis  of  order  two.  Again,  write  down  the  corresponding 
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matrices,  but  this  time  as  both  row  (one-by-two)  and  column 
(two-by-one)  matrices,  as  shown  below. 


(+1  +1)  (-1  +1) 


Multiply  together  the  two  numbers  located  on  the  outer 
row-column  intersections,  and  collect  the  products  in  a  two- 
by-two  groupings,  as  shown.  Note  that  these  groupings,  called 
Kronecker  products,  are  exact  numerical  analogs  of  the  mask 
pixel  structures  in  Figure  1.  The  new  multiplexed  imaging 
technique  will  implement  an  optical  scheme  to  incorporate  this 
property,  thereby  greatly  reducing  the  number  and  complexity 
of  the  encoding  mask  families  required  [Ref.  1].  Figure  2 
illustrates  a  complete  64  pixel  "checkerboard”  mask  basis, 
synthesized  from  two  one-dimensional  bases  of  order  eight. 

D.  RESTATEMENT  OF  THESIS  GOALS  IN  WALSH  SEQUENCY-THEORETICAL 

TERMS 

As  a  consequence  of  the  foregoing  discussion,  the  two 
goals  of  this  thesis  research  expressed  in  section  I.C.  may  be 
rephrased  in  sequency-theoretical  terms:  (1)  since  the 
prototype  system  will  generate  elaborate  two-dimensional  image 
encoding  masks  via  the  Kronecker  product,  the  system  designers 
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Figure  2  Two-dimensional  Seguency-Ordered  Walsh  Masks 

for  64  Pixels 
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need  to  know  how  these  masks  will  diffract  long  wavelength 
radiation,  as  well  as  encode  it.  Diffraction  effects  could 
cause  severe  beam  spreading  problems,  resulting  in  substantial 
loss  of  optical  efficiency,  particularly  at  these  longer 
wavelengths  and/or  when  using  encoding  masks  of  high  spatial 
sequency.  Therefore,  the  first  goal  of  this  research  has  been 
to  predict  these  effects  quantitatively.  (2)  The  second  goal 
is  stated  more  succinctly:  the  system  designers  need  to  know 
how  many  spatial  sequencies  are  sufficient  to  encode  images 
which  remain  recognizable  once  decoded.  It  would  be  both 
expensive  and  more  optically  complex  to  build  a  prototype 
instriiment  with  a  far  larger  Walsh  basis  set  (i.e.,  containing 
bases  of  higher  order)  than  actually  needed. 
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III.  THE  DIFFRACTION  PROBLEM 


As  noted  in  the  Chapters  I  and  II,  one  of  the  basic  tasks 
of  this  thesis  research  project  was  to  simulate  the  behavior 
of  the  proposed  multiplexed  imaging  scheme  at  long 
wavelengths.  Therefore,  this  first  phase  of  the  research  is 
called  the  diffraction  problem.  It  seeks  to  quantify  the 
diffraction  behavior  of  the  encoding  masks  by  means  of 
computer  simulation. 

A.  BASIC  SCALAR  DIFFRACTION  THEORY 

No  discussion  of  an  imaging  system  is  complete  without 
considering  the  effects  of  interference  and  diffraction  of  the 
incident  wavefronts.  The  importance  of  these  effects  cannot 
be  understated. 

"Diffraction  effects  are  accordingly  of  great 
significance  in  the  detailed  understanding  of  devices 
containing  lenses,  stops,  source  slits,  mirrors,  and  so 
on.  If  all  defects  in  a  lens  system  were  removed,  the 
ultimate  sharpness  of  the  image  would  be  limited  by 
diffraction."  [Ref.  3] 

The  essential  feature  of  diffraction  is  that  an  obstruction 
deviates  the  rectilinear  propagation  of  a  light  wave.  The 
geometry  of  diffraction  through  a  rectangular  aperture  is 
summarized  by  Hecht  [Ref.  3]  and  depicted  in  Figure  3. 
References  to  the  aperture  plane  are  in  lower  case  letters; 
those  to  the  image  plane  are  in  capitals. 
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Figure  3  Diffraction  Through  a  Rectangular  Aperture 
As  predicted  by  the  Huygens-Presnel  principle,  a 
differential  element  of  surface  area  within  the  aperture,  dS 
(equivalent  to  dydz) ,  may  be  considered  as  comprised  of  an 
infinite  number  of  secondary  wave  sources.  However,  for  a 
source  of  monochromatic  light,  dS  is  much  smaller  than  the 
wavelength  A,  and  all  of  dS's  contributions  to  the 
electromagnetic  disturbance  at  point  P  in  the  image  plane 
arrive  in  phase  and  interfere  constructively.  Assuming 
uniform,  normally  incident  illumination  over  the  entire 
aperture,  it  follows  that 
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E(ky,k^)  «  J  J  A(y,z)  dydz 


(5) 


where  the  aperture  function  is 


A(v  z'\  -  I  ^  (y»2)  -i®  transparent  regions 

"  I  0  for  (y,z)  in  opaque  regions 

(6) 

The  wave  number  k,  is  2n/k,  The  spatial  frequencies  in  the  y 
and  z  directions  are,  respectively. 


ky  -  xY/R  -  Ksin0  , 

(7) 


and 

ICj  -  x2/R  •  KsinO  . 

(8) 

E  represents  the  diffracted  electric  field  distribution. 
Figure  4  illustrates  the  angular  relationship  of  0  and  e 
between  the  aperture  and  image  planes. 

Equation  (5)  indicates  that  the  electromagnetic  field 
dist-' ■  bution  at  the  image  plane  is  proportional  to  the  two- 
dimensional  Fourier  transform  of  the  aperture  function. 
Determining  the  diffraction  behavior  of  the  photon  flux 
through  the  Walsh  masks  requires  a  three-step  calculation. 
First,  simulate  the  aperture  function  for  each  mask  of 
interest.  Second,  calculate  the  two-dimensional  Fourier 
transform  of  that  mask.  Third,  calculate  the  actual 
distribution  of  the  resulting  irradiance  as  a  function  of 
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X 


Figure  4  Diffraction  Geometry 

spatial  frequencies  ky  and  kj.  The  third  step  is  simple;  the 
irradiance  I  is  proportional  to  the  time  average  of  |E|^ 


[Ref.  3],  such  that 


I  •ci  EE*} 


(9) 


B.  OENERATIOM  07  mLBE  MASK  APERTURE  FUNCTIONS 

The  orthogonal  Walsh  aperture  functions  were  simulated  by 
means  of  a  straightforward  algorithm,  which  is  easily 
described  via  an  analogy  with  the  familiar  topic  of  Fourier 
transforms.  To  recover  one  of  the  orthogonal  trigonometric 
basis  functions  upon  which  Fourier  analysis  is  built,  one 
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could,  in  principle,  calculate  the  Fourier  transform  of  a 
Dirac  delta  function  pair  [Ref.  5],  such  that 

r  [5(/-/o)  +  5(/+fo)  df  - 

-  2  cos(2irfot} 

(10) 

vniile  there  are  more  efficient  methods  of  calculating  a  simple 
cosine  function,  the  transform  approach  is  conceptually  soiuid. 
Similarly,  if  one  has  at  their  disposal  an  efficient  Walsh 
transform  routine,  it  becomes  a  trivial  matter  to  generate 
Walsh  functions  of  different  seguencies:  simply  Walsh 
transform  a  nvimerical  ecpiivalent  of  a  delta  function  at  the 
desired  seguency. 

The  mathematical  construction  of  the  Walsh  masks  as 
outlined  above  was  accomplished  by  Davis  using  a  fast  Walsh 
transform  algorithm  which  he  developed  from  the  generalized 
transform  algorithm  theory  presented  in  Elliot  and  Rao  [Ref. 
7].  This  program  was  originally  written  in  Microsoft  C. 
Davis  translated  the  code  into  FORTRAN  for  incorporation  into 
the  mask  diffraction  program  WMASK.for  (listed  in  Appendix  A) 
as  a  subroutine  called  FWT.  The  necessary  parameters  that 
govern  the  mask  generation  and  ultimate  diffraction  pattern 
are  the  desired  spatial  seguency,  the  size  of  the  basis  set  of 
Walsh  functions,  and  the  overall  dimensions  of  the  two- 
dimensional  mask  matrix  which  will  be  produced.  For  this 
stage  of  the  project,  it  was  decided  that  the  number  of 
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sequencies  in  the  basis  set  would  remain  fixed  at  64, 
corresponding  to  an  eight-by-eight  array  of  square  pixels,  as 
shown  in  Figure  2.  The  interactive  portion  of  program 
WMASK.for  requests  specification  of  the  desired  mask  sequency 
in  the  range  zero  to  63. 

The  FWT  subroutine  operates  on  a  delta  function  and 
returns  a  one-dimensional  array  of  64  elements  corresponding 
to  the  discretely  sampled  Walsh  function  with  the  specified 
sequency  number.  These  64  elements  have  assigned  values  of  +1 
or  -1,  representing  the  transparent  and  opaque  regions, 
respectively,  of  the  Walsh  masks.  Thb  program  adds  the  value 
one  to  each  returned  element,  and  then  divides  each  element  by 
two.  This  normalizes  the  sequency  array  to  values  of  either 
one  (for  a  transparent  pixel)  or  zero  (for  an  opaque  pixel) . 
As  a  visual  check,  whose  importance  was  later  realized,  the 
program  writes  the  64  elements  to  the  computer  screen  in  the 
form  of  a  eight -by-eight  matrix.  Trial  and  error  simulations 
and  numerical  experiments  revealed  that  this  basic  eight-by¬ 
eight  matrix  format  was  undersampled  for  subsequent  Fourier 
transformation  [Ref.  5],  so  the  matrix  structure  was  expanded 
by  a  factor  of  two  in  each  dimension.  For  clarification,  the 
following  illustrates  the  process  in  generating  the  Walsh  mask 
for  sequency  28  of  63  (the  16-square  checkerboard  in  Figure 
2)  .  The  64  element  array  returned  by  the  subroutine  FWT  is 
ordered  into  an  eight-by-eight  matrix  in  the  form 
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11001100 

11001100 

00110011 

00110011 

11001100 

11001100 

00110011 

00110011 

After  the  bi-directional  expansion  by  a  factor  of  two,  the 
eight-by-eight  matrix  becomes  a  16  by  16  matrix  appearing  on 
screen  as 

1111000011110000 

1111000011110000 

1111000011110000 

1111000011110000 

0000111100001111 

0000111100001111 

0000111100001111 

0000111100001111 

1111000011110000 

1111000011110000 

IIIIOOOOIIIIOOOO 

1111000011110000 

0000111100001111 

0000111100001111 

0000111100001111 

0000111100001111 


The  two-dimensional  Fourier  transform  subroutine  used  in 
the  program  WMASK.for  mandates  that  the  input  matrix  contain 
complex  numbers.  The  required  format  calls  for  the  real  part 
of  the  number  followed  by  the  imaginary  part,  in  column  order. 
However,  the  elements  which  comprise  the  Walsh  mask  are  real 
numbers.  For  this  reason,  the  numerical  model  of  the  Walsh 
mask  as  shown  above  must  be  modified  by  injecting  a  row  of 
zeros,  representing  the  imaginary  component  of  each  element. 


20 


between  each  row.  This  transforms  the  mask  to  the  proper 
complex  representation.  The  resulting  matrix  is  32  by  16. 

Finally,  the  mask  was  centered  on  a  64  by  32  zero  matrix 
in  complex  form,  which  constitutes  a  border  of  zeros,  or 
zeropad,  of  a  factor  of  two  surrounding  the  mask.  This 
dimension  doubling  was  again  reached  by  trial  and  error.  Its 
purpose  was  to  interpolate  the  resulting  spectrum  via 
oversampling  [Ref.  5].  By  increasing  the  overall  dimensions 
of  the  input  matrix  by  a  factor  of  two  (reaching  final 
dimensions  of  32  by  32) ,  the  inherent  spectral  sampling  rate 
is  doubled,  and  consequently,  the  resulting  spatial 
frequencies  in  each  dimension  of  the  diffraction  pattern  are 
then  more  densely  spaced  by  a  factor  of  two. 

As  outlined  in  Press,  et  al.,  [Ref.  6],  the  complex  data 
must  first  be  properly  ordered  prior  to  performing  the  two- 
dimensional  Fourier  transform.  This  ordering  is  required  to 
reduce  the  number  of  steps  necessary  to  complete  the 
transformation.  As  is  typical  of  many  one-dimensional  FFT 
programs,  the  input  data  must  be  reordered  in  a  wrap-around 
fashion  by  placing  the  second  half  of  the  data  in  front  of  the 
first  half  prior  to  conducting  the  transform.  This  quadrant 
reversal  is  also  discussed  in  Brigham  [Ref.  5].  In  two 
dimensions,  with  the  data  in  matrix  form,  this  wrap-around 
procedure  must  be  performed  twice,  once  in  each  spatial 
direction.  Finally,  the  matrix  is  passed  to  the  two- 
dimensional  FFT  subroutine  in  normal  FORTRAN  order.  This 
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ordering  requires  partitioning  the  matrix  into  a  one¬ 
dimensional  array  composed  of  the  sequence  of  the  vertical 
columns  of  the  input  matrix.  Such  an  ordering  structure  is 
well  represented  in  Figure  12.11.1  of  Press,  et  al.,  [Ref.  6]. 
The  FFT  subroutine  also  returns  the  transformed  data  in  this 
order,  which  must  be  sorted  following  the  reverse  of  the  input 
procedure,  to  recover  the  two-dimensional  FFT  of  the  input 
matrix. 

The  following  illustrates  the  steps  involved  in  the  data 
manipulation  of  a  simple  two-by-two  input  matrix  enroute  to 
the  FFT  subroutine.  The  numbers  represent  the  matrix 
elements.  Each  step  was  originally  developed  as  a  separate 
program  to  verify  its  function,  and  then  inserted  into  the 
master  program  WMASK.for  as  subroutines. 


1. 

2. 

3. 


1  2  Input  Data  Matrix 

3  4 

1  2  Inject  Zeros  to  Transform  to  Complex  Data 

0  0 
3  4 

0  0 

0000  Center  Complex  Matrix  on  Zeropad 
0  0  0  0 
0  12  0 
0  0  0  0 
0  3  4  0 

0  0  0  0 
0  0  0  0 
0  0  0  0 
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4.  0000  Horizontal  Wrap-Around 

0  0  0  0 

2  0  0  1 
0  0  0  0 
4  0  0  3 

0  0  0  0 
0  0  0  0 
0  0  0  0 

5.  4003  Vertical  Wrap-Around 

0  0  0  0 

0  0  0  0 
0  0  0  0 
0  0  0  0 
0  0  0  0 
2  0  0  1 
0  0  0  0 

6.  4  0  0  3  Form  1-D  Array  to  Subroutine 

0  0  0  0 

0  0  0  0 

0  0  0  0 

0  0  0  0 

0  0  0  0 

2  0  0  1 

0  0  0  0 


In  returning  the  transformed  data,  the  steps  are  taken  in 
reverse . 


C.  THE  FAST  FOURIER  TRANSFORM 

The  fast  Fourier  transform  (FFT)  has  played  a  significant 
role  in  reducing  the  amount  of  computing  time  required  for 
Fourier  analysis.  In  two  or  more  dimensions,  as  is  required 
by  this  and  other  imaging  projects,  the  nvimber  of  calculations 
becomes  so  large  that,  without  the  FFT,  it  could  not  be  done 
in  a  useful  amount  of  time.  The  key  feature  which  makes  the 
FFT  work  so  efficiently  is  an  algorithm  to  reduce  redundant 
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arithmetic  operations.  Brigham  [Ref.  5]  outlines  the  theory 
of  the  FFT  and  describes  its  implementation. 

The  two-dimensional  FFT  used  in  this  project  is  an 
adaptation  of  one  taken  from  Press,  et  al.,  [Ref.  6],  which 
was  developed  by  N.  M.  Brenner  of  Lincoln  Laboratories.  The 
program  provided  in  the  reference  was  written  in  FORTRAN. 
Consequently,  all  programming  for  this  part  of  the  thesis  was 
done  in  this  language.  Although  several  personal  computers 
were  used  throughout  the  programming  phase,  the  majority  of 
the  work  was  performed  on  a  16-bit  Imperial  80286  machine, 
using  the  Microsoft  FORTRAN  4.1.  compiler.  To  minimize 
computing  time,  an  executable  file  of  the  successfully  tested 
program  was  installed  on  an  IBM  Model  70  PS/2  80386  personal 
computer  with  an  80387  coprocessor. 

0.  CONVERTING  THE  TRANSFORMED  MATRIX  INTO  DIFFRACTION 

PATTERNS 

The  array  returned  from  the  FFT  subroutine  is  reordered 
back  into  complex  matrix  form,  but  does  not  yet  represent  the 
diffraction  pattern  of  the  input  Walsh  mask.  In  reference  to 
equation  (5) ,  the  two-dimensional  Fourier  transform  of  the 
aperture  function  is  an  expression  for  the  electric  field 
distribution  E(ky,k2)  at  the  Fourier  plane.  However,  the 
desired  quantity  is  the  irradiance  I,  as  a  function  of  the 
spatial  frequencies. 
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The  algorithm  for  implementing  equation  (9)  is  contained 
within  the  subroutine  INTENSITY  of  program  WHASK.for.  This 
subroutine  converts  the  returned  complex  field  distribution 
matrix  into  a  real  matrix  where  each  element  represents  the 
intensity  at  a  point  in  the  Fourier  diffraction  pattern.  In 
doing  so,  the  matrix  is  reduced  from  dimensions  of  64  by  32, 
to  a  32  by  32  square.  The  matrix  then  represents  the  field 
intensity  as  a  function  of  spatial  frequency,  which  is  the 
desired  form  for  the  diffraction  pattern  of  the  Walsh  mask 
aperture  function. 

A  consequence  of  performing  a  two-dimensional  fast  Fourier 
transform  on  a  square  matrix  is  that,  in  the  final  form,  the 
first  row  and  the  first  coliunn  contain  all  zeros.  This  result 
stems  from  the  vertical  and  horizontal  quadrant  reversal 
required  in  the  matrix  manipulation  prior  to  transformation. 
As  this  first  row  and  column  have  no  real  meaning,  they  are 
eliminated  from  the  intensity  matrix.  Therefore,  this  last 
operation  leaves  the  matrix  in  the  form  of  a  31  by  31  square. 

E.  PLOTTING  THE  DIFFRACTION  PATTERNS 

It  was  necessary  to  locate  a  software  package  which  could 
produce  a  three-dimensional  perspective  plot  of  the  intensity 
matrix  as  a  function  of  spatial  frequency,  such  a  program  is 
SURFER,  which  is  distributed  by  Golden  Software,  Inc.  This 
program  was  used  to  produce  all  of  the  diffraction  pattern 
plots  which  appear  in  Appendix  B. 
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The  format  required  by  SURFER  demands  that  the  data  to  be 
plotted  be  accessible  in  three-dimensional  Cartesian  (x,y,z) 
coordinates.  At  this  stage  in  the  programming,  however,  the 
intensity  was  represented  as  a  31  by  31  square  matrix.  A 
subroutine  called  MAT2XyZ  was  developed  to  convert  the  matrix 
into  the  required  format.  In  this  form,  each  converted  point 
was  assigned  a  set  of  (x,y)  coordinates  to  locate  the  point  in 
the  intensity  matrix.  The  (z)  coordinate  indicated  the 
intensity  at  that  point.  Such  a  process  generates  961  data 
points  from  the  31  by  31  intensity  matrix,  and  writes  them  to 
a  disk  file  for  later  use  by  the  plotting  routine. 

SURFER  employs  a  two-step  process  to  create  three- 
dimensional  plots  of  input  data.  The  first  step  involves 
forming  a  regularly  spaced  grid  from  irregularly  spaced  data. 
The  data,  in  the  case  of  this  project,  are  the  (x,y,z) 
elements  which  are  generated  by  the  master  program  as 
described  above.  These  data  points  are  imported  from  the  disk 
file  into  SURFER  by  specification  during  the  gridding  phase. 
In  the  language  of  the  software,  this  process  is  simply  called 
"gridding”,  and  is  invoked  in  the  menu  heading  GRID.  The 
second  step  in  the  plotting  process  calls  for  loading  the 
gridded  data  into  the  subprogram  called  SURF.  SURF  is  a  menu- 
driven,  interactive  graphics  routine  which  produces  the 
perspective  plots  constituting  the  diffraction  patterns.  Its 
output  is  a  plot  file  which  contains  the  SURFER-generated 
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surface  plot  data.  This  plot  file  may  be  viewed  on  screen  or 
sent  directly  to  a  specified  plotter. 

As  noted,  SURFER  is  a  menu-driven  program.  The  following 
parameters  were  used  during  the  interactive  set-up  to  create 
the  diffraction  pattern  plots  which  appear  in  Appendix  B.  The 
parameters  are  listed  in  the  order  in  which  they  were  input, 
and  are  the  result  of  much  trial  and  error  to  find  the  optimal 
settings.  Many  of  the  default  parameters  were  satisfactory 
for  making  the  diffraction  pattern  plots;  only  the  parameters 
which  differ  from  the  defaults  are  listed  below. 


At  the  Main  Menu:  GRID 

Filespec  to  Pass:  XYZ.GRD 

At  Menu :  Modify 

Select:  Smooth 

Select:  Spline  (expansion  factor  2) 

Select:  Begin 


This  procedure  creates  the  output  grid  file  OUT.GRD 


At  the  Main  Menu: 

At  Menu: 

Input  grid  file: 

At  Menu: 

Projection: 

At  Menu: 

Plot  constant  (x,y,2)  : 
At  Menu: 

Plot  base: 

At  Menu: 

Plot  Title 

Title  symbol  set: 

Title  position: 

Title  angle  in  degrees: 
Title  character  height: 
Title  color: 

Plot  orientation  legend: 
Legend  position: 


SURF 

Input 

OUT.GRD 

View 

Perspectivo 

LineTyp 

Base 

Mo 

Title 

Diffraction  Pattern 
Tor  Walsh  Mask  #28 
DEFAULT. 8 YN 
(0,-1) 

0 

0.15 

1 

NO 

Automatic 
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At  Menu: 

Plot  axes: 

Axis  color: 

Axis  symbol  set: 

Axis  plane: 

Tic  distance: 

Labeled  tic  frequency 
Label  format: 

Number  of  decimal  digits 
Label  character  height: 
Tic  label  angle: 

Axis  title: 

Title  character  height: 
Title  to  label  distance: 

No  changes  are  required  in  menu 

Post. 

At  Menu: 

Name  of  plot  file: 

Scale  factor: 

Page  position: 

Plot  file  format: 

File  write  mode: 

Number  of  decimal  digits: 
Send  to  output  device: 

No  changes  are  required  in  menu 


Axes 

Z,Y  (each  separately) 

1 

DEFAULT. 8 YM 
Automatic 
6.28 
2 

Fixed 

2 

0.1 

0 

Spatial  Frequency  X  sin  A/B 

0.13 

0.14 

headings  Size,  ZYLine,  or 


Output 

OUT.PLT 

1 

1/1 

Binary 

Overwrite 

3 

Mo 

heading  Equip. 


This  procedure  generates  the  plot  file  OUT.PLT,  which  is 
called  for  in  the  final  phase  of  SURFER.  This  file  is  sent  to 
the  HP7550A  plotter. 


At  Main  Menu: 

Name  of  plot  file: 

Do  you  wish  to  shift 
Do  you  wish  to  scale 
Scale  factor  in  X 
Scale  factor  in  Y 


PLOT 
OUT . PLT 

the  entire  plot?  No 
the  entire  plot?  Yes 
direction:  0.6 

direction :  0.6 
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F.  ANALYSIS  OF  DIFFRACTION  PATTERNS 


Having  established  a  working  diffraction  program  and 
stored  the  listed  plotting  parameters  into  SURFER,  the 
remainder  of  this  phase  of  the  thesis  research  was  devoted  to 
generating  the  Walsh  masks  and  to  examining  their  diffraction 
pattern  plots  to  reach  conclusions  about  their  diffraction 
behavior.  Analysis  of  the  first  few  diffraction  pattern  plots 
immediately  raised  guestions  about  what  they  represented;  the 
results  did  not  seem  to  conform  to  expectations.  The  computer 
code  was  reviewed  to  determine  where,  if  any,  mistakes  were 
present.  "What  should  the  diffraction  patterns  of  the  Walsh 
masks  should  look  like?"  became  the  question  of  interest. 
Consequently,  the  project  focus  shifted  to  resolve  this 
matter . 

G.  APPARENT  PROBLEMS  WITH  THE  DIFFRACTION  PATTERNS 

The  first  few  diffraction  patterns  plotted  were  very 
encouraging.  In  reference  again  to  the  Walsh  masks  in  Figure 
2,  sequency  zero  appears  in  the  upper  left  hand  corner  as  the 
solid  black  square.  When  centered  on  the  zeropad,  the  black 
square  is  seen  as  a  square  aperture  (considering  the  black 
portion  as  transparent,  the  zeropad  as  opaque) .  The 
diffraction  pattern  of  a  square  aperture  is  well  known  and 
very  identifiable.  It  appears  as  Figure  10.25(b)  in  Hecht 
[Ref.  3].  It  was  with  much  relief  that  the  plot  generated  by 
SURFER  of  Walsh  sequency  zero  (Figure  B-1)  was  identical  to 
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this  figure  within  the  scaling  limits.  Continuing  across  the 
first  row  of  masks «  with  seguencies  one,  two,  etc.,  the 
diffraction  patterns  bore  the  features  which  would  have  been 
expected  for  an  array  of  vertical  slits.  Despite  the  apparent 
success,  the  diffraction  pattern  produced  by  seguency  63  (the 
64-sguare  checkerboard  in  Figure  2)  did  not  conform  at  all  to 
expectations. 

This  checkerboard  pattern  has  the  finest  detail  of  all  of 
the  masks  in  the  basis  set  of  Figure  2.  Conseguently,  its 
diffraction  pattern  would  be  expected  to  be  spread  out  much 
more  than  the  others.  Also,  its  symmetry  about  the  main 
diagonal  predicts  a  likewise  symmetric  diffraction  pattern. 
The  plot  generated  from  seguency  63  was  neither  well  dispersed 
nor  symmetric.  Only  after  installing  the  routine  to  display 
the  mask  on  the  monitor  did  the  explanation  for  this  become 
apparent. 

1.  False  Assumptions  in  the  Sequenoy  Ordering 

From  the  outset,  it  was  assumed  that  the  Walsh  masks  of 
Figure  2  were  seguentially  ordered  left-to-right  in  a  raster 
fashion,  beginning  at  the  upper  left  hand  corner.  Certainly, 
the  black  sguare  represents  seguency  zero.  In  this  system 
then,  the  masks  of  this  figure  appear  "numbered"  in  seguency 
as 
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01234567 
8  9  10  11  12  13  14  15 
16  17  18  19  20  21  22  23 
24  25  26  27  28  29  30  31 
32  33  34  35  36  37  38  39 
40  41  42  43  44  45  46  47 
48  49  50  51  52  53  54  55 
56  57  58  59  60  61  62  63 

In  such  a  system,  the  first  row  appears  in  order,  and  the 
checkerboard  is  located  in  position  63. 

The  revelation  caused  by  the  routine  to  display  the  masks 
indicated  an  entirely  different  ordering.  By  calling  for  the 
masks  in  ascending  order  (from  zero),  it  became  evident  that 
instead  of  a  raster  structure,  the  subroutine  FWT  returned  the 
data  in  a  serpentine  fashion.  Regardless  of  this 
understanding,  the  first  row  appears  identically  in  each 
scheme.  This  explains  why  the  scpiare  aperture  pattern  could 
be  generated  without  full  comprehension  of  the  sequency 
ordering.  However,  the  mistakenly  assumed  raster  ordering 
places  the  checkerboard  in  a  different  location  as  opposed  to 
where  this  mask  properly  appears  when  located  with  a 
serpentine  ordering.  Therefore,  rather  than  generating  the 
diffraction  pattern  of  the  checkerboard  by  specifying  sequency 
63,  the  diffraction  pattern  of  the  mask  located  in  position  56 
was  being  produced,  which  is  predictably  neither  well 
dispersed  nor  symmetric. 
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2.  Tb«  Proper  Ordering  of  the  Welsh  Masks 

with  the  preceding  discussion  in  mind,  the  proper 
"numbering"  of  the  Walsh  masks  as  generated  by  the  diffraction 
pattern  program  and  appearing  in  Figure  2  is 

01234567 

15  14  13  12  11  10  9  8 

16  17  18  19  20  21  22  23 

31  30  29  28  27  26  25  24 

32  33  34  35  36  37  38  39 

47  46  45  44  43  42  41  40 

48  49  50  51  52  53  54  55 

63  62  61  60  59  58  57  56 

Throughout  the  remainder  of  this  thesis,  the  Walsh  masks  and 
their  subsequent  diffraction  patterns  will  be  referred  to 
using  this  system,  whereby  the  assigned  number  indicates  both 
the  sequency  of  the  mask  and  its  location  in  the  above  table. 

H.  SYMMETRIES  IN  THE  WELSH  MESKS 

With  a  complete  understanding  of  the  ordering  of  the 
masks,  it  was  quickly  determined  that  the  progreun  had  been 
working  all  along,  and  production  of  the  diffraction  pattern 
plots  resumed  in  earnest.  However,  not  all  of  the  masks  were 
used.  Collectively,  an  examination  of  the  masks  in  Figure  2 
reveals  that  masks  reflected  across  the  main  diagonal  are 
identical  except  for  a  90  degree  rotation.  This  symmetry 
allows  one  to  draw  conclusions  about  the  complete  set  of  masks 
by  examining  only  the  upper  (or  lower)  triangle  of  the  masks 
in  this  figure.  Such  a  symmetry  reduces  the  number  of 
required  analyses  and  plots  from  64  down  to  36,  as  the 
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dll  fraction  pattern  of  a  rotated  mask  is  the  same  as  its  non- 
rotated  counterpart;  the  only  difference  being  that  the 
pattern  is  also  rotated  by  the  same  amount  as  the  mask.  The 
diffraction  patterns  which  comprise  Appendix  B,  then,  are 
taken  from  the  36  masks  which  appear  on  and  above  the  main 
diagonal  of  Figure  2. 

I.  A  TYPICAL  DIFFRACTION  PATTERN  PLOT 

Figure  5  shows  the  diffraction  pattern  for  Walsh  mask  28. 
This  mask  is  the  16-square  checkerboard.  This  typical  pattern 
illustrates  the  type  of  information  which  is  analyzed  when 
examining  these  plots.  Angles  A  and  B  correspond  to  angles  <p 
and  e  shown  in  Figure  4;  they  relate  to,  respectively,  the 
angular  spread  due  to  diffraction  in  the  Z  and  Y  directions  of 
the  Fourier  plane.  The  axes  are  labeled  in  normalized  units 
of  spatial  frequency  corresponding  to  ky  and  kj.  Each 
direction  is  labeled  out  to  values  of  ±  2n.  The  multiples  of 
n  are  used  to  keep  the  dimensions  consistent  in  terms  of 
spatial  frequency.  These  units  of  spatial  frequency  are  the 
reciprocal  of  the  Walsh  mask  pixel  size. 

The  irradiance  (vertical)  axis  is  neither  labeled  nor 
scaled.  This  fact  is  due  to  the  nature  of  the  Walsh  masks. 
All  masks  (except  for  the  square  aperture  of  sequency  zero) 
have  equal  transparent  and  opaque  areas.  Therefore,  the 
diffraction  patterns  of  all  of  the  masks  will  transmit  the 
same  amount  of  electromagnetic  energy,  albeit  in  different 
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distributions.  Hence,  the  relative  irradiance  scaling  is  the 
same  for  all  masks,  and  absolute  scaling  is  unnecessary. 

J.  CONCLUSIONS  CONCERNING  THE  DIFFRACTION  LIMIT 

The  quantity  which  concerns  us  is  the  spreading  of  the 
irradiance  in  each  diffraction  pattern.  Although  all  of  the 
patterns  show  a  characteristic  central  peak  containing  the 
majority  of  the  irradiance,  some  display  a  significant  amount 
of  energy  spilling  out  into  the  periphery  of  the  pattern.  For 
reference  purposes.  Appendix  B  constitutes  a  catalog  of  all  36 
different  patterns  generated  by  the  Walsh  masks  in  the  basis 
set  of  64.  They  may  be  examined  individually  to  determine  the 
limit  in  terms  of  the  spatial  frequency  necessary  to  retain 
the  diffracted  irradiance;  or,  perhaps  more  expeditiously, 
only  the  extreme  cases  need  be  considered.  The  result  of  such 
a  process  constitutes  the  answer  to  the  research  question  at 
hand. 

In  reference  to  the  Walsh  masks  in  Figure  2,  the  finest 
detail  occurs  in  mask  56,  the  dense  checkerboard  pattern. 
Certainly  its  diffraction  pattern  must  be  considered  as  a 
limiting  case.  Nearly  all  of  the  masks  in  the  first  row  show 
diffraction  patterns  with  wide  angular  dispersion  as  well, 
particularly  masks  1,  2,  and  3.  Such  patterns  define  the 
extremes  of  diffraction  spreading  for  this  basis  set  of  Walsh 
masks.  It  is  important  to  ensure  that  the  limit  chosen 
encloses  nearly  all  of  the  significant  irradiance  of  even  the 
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most  severe  case  of  diffraction.  After  careful  consideration, 
it  was  determined  that  this  worst  case  occurs  for  Walsh 
seguency  three  (Figure  B-4) . 

K.  FORMULATING  THE  DIFFRACTION  LIMIT 

Having  made  the  above  determination,  it  became  a  matter  of 
simply  reading  from  the  plot  the  minimum  spatial  frequency 
required  to  encompass  the  majority  of  the  diffracted 
irradiance.  It  was  concluded  that  this  minimum  spatial 
frequency  was  6  reciprocal  length  units,  as  defined  above.  To 
put  this  result  in  a  form  useful  in  future  instrumental 
design,  it  may  be  stated  that  the  following  relation  must  be 
satisfied 

1^  sindl  s  6  . 

(11) 

The  angle  e  is  used  here  since,  once  a  mask  is  "rotated",  in 
that  its  counterpart  across  the  main  diagonal  is  used,  the 
angle  A  becomes  equivalent  to  the  angle  B.  For  simplicity, 
the  diffraction  angle  henceforth  shall  be  called  0.  Figure  6 
illustrates  the  simplified  relationship  between  the  mask,  the 
optics,  and  the  detector. 
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Mask 


Optics 


Detector 


Figure  6  Mask,  Optics,  and  Detector  Orientation 
The  formulation  of  equation  (11)  constitutes  the  solution 


to  the  first  thesis  problem. 
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IV.  RESOLUTION  CRITERIA  FOR  IMAGE  PATTERN  RECOGNITION 


A.  OVERVIEW 

The  second  phase  of  this  research  project  involved,  by 
means  of  computer  simulation,  a  determination  of  what  minimal 
spatial  resolutions  are  needed  for  pattern  recognition  in 
images  of  commonplace  objects.  Establishment  of  these 
criteria  are  of  great  importance  in  the  overall  multiplexed 
imaging  research  project,  for  they  dictate  the  minimal  number 
of  picture  elements  (pixels)  needed  by  a  practical  remote 
sensing  instrument.  This  minimal  number  of  pixels,  in  t\irn, 
specifies  the  number  of  Walsh  basis  functions  and  encoding 
masks  required  to  accomplish  this  task. 

A  working  prototype  system  does  not  yet  exist;  therefore 
this  problem  was  addressed  in  reverse.  Digitized  images  were 
transformed  into  a  sets  of  sequency  spectra .  These  spectra 
were  sequency-truncated  (low>pass  filtered)  to  mimic  the 
effects  of  limited  resolution  in  an  actual  instrument.  The 
filtered  sequency  spectra  were  then  re-transformed  into  the 
image  domain  for  analysis. 

The  procedure  may  be  further  clarified  by  recalling  that 
an  image's  higher  spatial  frequencies  are  responsible  for 
finer  image  detail.  In  this  research,  spatial  sequencies  play 
the  role  normally  reserved  for  spatial  frequencies  with  regard 
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to  traditional  image  analysis  approaches.  Elimination  of 
higher  seguencies  degrades  image  resolution  in  the  same  sense 
that  the  fine  detail  in  an  observed  field  is  degraded  by  the 
loss  of  components  of  higher-order  spatial  frequency. 

B.  THE  WALSH  FUNCTIONS  AND  IMAGE  PROCESSING 

Previous  sections  have  discussed  the  Walsh  functions  with 
regard  to  seguency  theory  and  the  proposed  encoding  technique. 
Beauchamp  [Ref.  4]  goes  much  further  in  detailing  several 
other  applications  of  the  Walsh  functions  in  image  processing, 
including  pattern  recognition,  and  image  enhancement, 
transmission,  and  restoration.  However,  those  properties  of 
the  Walsh  functions  which  are  summarized  in  Chapter  II  of  this 
thesis  are  sufficient  for  a  description  of  the  procedures  used 
in  this  phase  of  the  research. 

C.  CREATING  DIGITIZED  IMAGES 

The  actual  order  of  tasks  for  determining  the  nximber  of 
seguencies  required  for  addressing  the  image  resolution  - 
pattern  recognition  problem  were  as  follows: 

•  Capture  a  series  of  images  on  a  videocamera 

•  Digitize  the  images  into  computer  disk  files 

•  Write  software  to  store,  access,  and  display  the  images 

•  Write  software  to  perform  the  seguency  low-pass  filtering 

•  Analyze  the  resulting  images  to  establish  appropriate 
resolution  criteria. 
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D.  CAPTURING  THE  IMAGES 


The  original  plan  was  to  record  still  images  and  to 
perform  a  digitizing  scan  of  each  using  a  LYNXX  solid-state 
CCD  slow-scanning  imaging  camera.  Prof.  D.L.  Walters  procured 
such  a  camera  system  along  with  the  necessary  imaging 
software,  from  SpectraSource  Instruments.  However,  the  camera 
system  was  modified  for  higher  frame  rates,  and  was  shipped 
back  to  the  manufacturer.  Hence,  because  the  system  was  not 
available  for  use  in  this  project,  an  alternative  was  sought. 

The  most  readily  available  device  for  image  recording  is 
an  ordinary  VHS  camcorder.  A  Panasonic  VHS-C  Palmcorder  was 
borrowed  from  a  friend  to  record  the  images  used  in  this 
project.  A  series  of  15  images  were  recorded  from  areas  of 
interest  around  the  Monterey  Peninsula,  with  the  understanding 
that  the  VHS  tape  would  be  played  back  in  a  freeze-frame  mode 
into  digitizing  equipment  operated  by  the  NPS  Physics 
Department.  The  output  from  this  equipment  would  then  be  a 
digitized  disk  file  of  each  image  in  simple  binary  format. 

E.  DIGITIZING  THE  IMAGES  FROM  VHS  TAPE 

Difficulties  arose  with  regard  to  gaining  access  to  the 
original  image  digitization  equipment,  due  to  the  author's 
lack  of  familiarization  with  it.  Therefore,  a  quest  to  locate 
alternate  means  to  complete  the  task  was  undertaken.  An  ideal 
system  would  grab  a  frame  from  the  VHS  playback,  scan  it, 
digitize  it,  and  then  write  the  image  to  an  binary  disk  file. 
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Several  systems  were  located  throughout  the  NPS  campus  which 
apparently  had  this  capability. 

1.  Apple  Macintosh  IZx  Programs 

Two  programs  in  use  at  the  NPS  Linear  Accelerator  facility 
were  offered  by  Prof.  X.K.  Maruyama,  and  each  was  thoroughly 
investigated  for  suitability  to  the  image  digitization  task. 
The  programs  are  IMAGE  1.29,  and  PIXELWORKS.  Both  programs 
run  on  an  Apple  Macintosh  Ilx  personal  computer,  and  have 
similar  capabilities  well  beyond  the  simple  needs  of  this 
project.  The  programs  are  able  to  manipulate  and  digitize  an 
image,  which  is  played  back  from  the  camcorder  into  the  video¬ 
in  port  of  the  Macintosh.  Difficulties  arose  at  this  stage  of 
the  process. 

The  programs  offer  a  wide  range  of  digitized  data  formats; 
in  particular,  TI.'F  .'and  simple  TIFF),  PICT,  MacPaint, 
Palette,  and  RAW,  are  available.  However;  binary  is  not. 
Substantial  effort  was  made  to  resolve  this  problem,  but 
despite  a  collaborative  effort,  no  simple  solution  was  found. 
As  an  alternative  approach,  Davis  suggested  that  if  the  format 
of  the  data  type  was  known,  he  could  then  write  a  program  to 
convert  the  data  into  binary  bytes.  Unfortunately,  the 
supporting  documentation  for  the  Macintosh  imaging  software 
was  not  detailed  enough  to  provide  this  information. 
Therefore,  further  investigation  into  this  translation  method 
were  halted,  and  alternate  digitization  means  were  sought. 
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2.  Other  Means  of  Digitization 

Earlier  failures  in  producing  a  digitized  image  disk  file 
in  binary  format  lead  to  an  expanded  search  for  an  alternate 
means  to  accomplish  this  task.  Several  Hewlett-Packard  page 
scanners  are  located  in  different  departments  at  NPS,  and  were 
considered  as  possible  candidates.  The  NPS  Computer  Center 
also  uses  HP  scanners,  as  well  as  an  IBM  page  scanning  system. 
The  IBM  system  appeared  most  promising,  but  it  proved  to 
demand  a  critical  investment  in  terms  of  training  time  and 
familiarization.  Furthermore,  it  was  perpetually 
oversubscribed.  Eventually,  these  page  scanners  were  dropped 
as  alternate  means  of  digitization. 

3 .  Successful  Digitisation 

Davis  contacted  Prof.  A.W.  Cooper  and  Mr.  W.J.  Lentz  of 
the  NPS  Physics  Department,  who  agreed  to  loan  some  of  their 
equipment  and  expertise  to  assist  the  author  with  the  image 
digitization.  Their  equipment  configuration  mandated  that  the 
original  VHS-C  camcorder  images  be  translated  onto  a  standard 
VHS  cassette  tape.  This  was  easily  done  using  the  playback 
feature  of  the  camcorder  in  conjunction  with  the  record  mode 
of  the  author's  VHS  recorder. 

The  digitizing  equipment  made  use  of  integrated  software 
coupled  through  an  IBM  personal  computer  to  link  a  playback 
VCR  and  a  display  monitor.  The  digitizing  program  was  DT  - 
IRIStutor  VOl.Ol,  which  is  marketed  by  DT,  Inc.  This  program 
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utilizes  an  IRIStutor  DT2861  frame-grabber  board  in 
conjunction  with  the  IBM  PC  to  reduce  a  captured  image  to  a 
512  by  512  pixel  digitized  image  file,  in  eight  bit  binary 
format.  The  images  from  the  VHS  tape  were  subsequently 
digitized  and  stored  on  5-1/4  inch  360k  byte  diskettes,  with 
one  image  per  diskette. 

F.  CONVERTING  AND  DISPLAYING  THE  IMAGES  ON  SCREEN 

The  files  generated  by  DT-IRIS  contained  a  header  block  of 
512  bytes  prior  to  the  actual  pixel-by-pixel  digitized 
representation  of  the  image.  It  was  necessary  to  strip  off 
this  information  prior  to  displaying  the  image  on  the  VGA 
monitor  of  the  project's  IBM  PS/2  model  70  PC.  In  order  to 
accomplish  this,  Davis  wrote  a  program  in  Microsoft  C5.1 
called  C512,  which  is  used  tc  convert  the  512  by  512  DT-IRIS 
images  into  simple  eight  bit  binary  format.  This  program  is 
included  as  Appendix  J. 

Once  the  image  file  had  been  converted  into  usable  form, 
it  was  necessary  to  write  another  routine  to  display  the  image 
on  the  monitor.  Again  programming  in  C,  Davis  developed  the 
program  D512,  which  accesses  the  binary  image  file,  and 
displays  it  in  16-level  greyscale  on  the  VGA  monitor.  In  this 
scale,  a  shade  of  zero  represents  black,  while  a  shade  of  15 
is  bright  white.  To  accommodate  the  line  restrictions  imposed 
by  the  monitor,  the  program  reduces  the  display  to  480  rows  by 
512  columns.  This  is  accomplished  by  stripping  off  the  top 
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16,  and  the  bottom  16,  horizontal  lines  of  the  image.  Thus, 
in  the  absence  of  any  seguency  filtering,  the  progreun  displays 
the  image  in  480  by  512  pixel  format.  The  program  D512  is 
included  as  Appendix  K. 

G.  FIVE  REPRESENTATIVE  IMAGES 

Originally,  a  total  of  15  scenes  were  recorded  using  the 
Panasonic  camcorder.  These  15  scenes  were  viewed  with  Davis 
to  select  those  which  contained  characteristic  features  and/or 
desirable  symmetries  which  might  prove  interesting  in  the 
context  of  the  low-pass  seguency  filtering  problem.  The 
following  images  were  chosen: 

•  U.S.  Coast  Guard  Cutter  Point  Barrow 

•  U.S.  Coast  Guard  Lighthouse  at  Point  Pinos 

•  A  friendly  German  shepherd 

•  A  house  with  a  picket  fence 

•  U.S.  Coast  Guard  emblem 

Despite  the  apparent  Coast  Guard  bias  in  subject  material, 
the  images  selected  represent  a  broad  range  of  spatial 
characteristics  such  as  fine  detail,  periodic  structure,  sharp 
contrast,  circular  symmetry,  and  so  forth.  The  original, 
unprocessed,  images  appear  in  the  Figures  7  through  11.  Note 
that  each  figure  is  captioned  with  a  noun  descriptor  followed 
by  a  nvimber.  This  naming  convention  has  been  adopted 
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throughout  the  remainder  of  the  research.  The  noun  name 
simply  describes  the  dominant  object  in  the  image;  the  number 
indicates  the  number  of  seguencies  used  in  each  direction  in 
the  image  basis  set.  Using  this  scheme,  the  first  image  is 
Cutter512,  indicating  that  it  is  the  Coast  Guard  cutter  with 
no  seguencies  eliminated  (i.e.,  512  by  512  pixels). 

This  labeling  system  is  used  throughout  the  project  and  is 
most  useful  when  studying  the  images  contained  within 
Appendices  C  -  I.  However,  there  is  one  departure  from  the 
system.  Images  labeled  with  a  “V"  or  an  "H”  postscript 
indicate  that  the  image  was  generated  using  only  a  one¬ 
dimensional  Walsh  transform,  with  that  transform  applied  in 
the  vertical  (V)  or  horizontal  (H)  direction.  These 
descriptors  appear  only  in  conjunction  with  the  seguency 
filtering  applied  to  the  image  of  the  cutter,  and  are 
contained  in  Appendices  C  and  D. 

All  of  the  images  which  appear  in  the  appendices  were 
photographed  directly  off  of  the  VGA  monitor.  For  this 
reason,  the  curvature  of  the  screen  is  noticeable.  The  small 
black  tab  which  seems  to  move  about  on  the  left  and  lower 
borders  is  a  marker  placed  by  the  author  to  identify  the 
developed  prints.  The  shaded  diagonal  strip  which  is  evident 
in  some  of  the  photographs,  such  as  in  Cutter2V  in  Figure  C-9, 
is  a  conseguence  of  the  monitor's  1/30  second  refresh  and 
flyback  time  when  photographed.  The  pictures  were  taken  with 
black  and  white  ISO  400  KODAK  Tri-X  pan  film,  at  an  aperture 
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setting  of  t4,  and  shutter  speed  of  1/15  second.  The  camera 
was  a  Canon  AE*-1. 

H.  PROCEDURES  BY  WHICH  IMAGE  RESOLUTION  WAS  MODIFIED 

The  FORTRAN  program  IMAGE. for,  written  by  the  author,  is 
the  primary  program  which  simulates  image  resolution 
degradation.  It  does  this  by  taking  an  original,  full 
resolution,  digitized  image  produced  by  program  C512,  and 
operates  on  it  to  generate  a  new  file  with  degraded 
resolution.  The  degradation  involves  three  steps.  First,  the 
program  calculates  the  Walsh  transform  of  the  full-resolution 
image,  to  produce  its  seguency  spectrum.  Second,  it  performs 
low-pass  seguency  filtering  on  the  image  spectrvun,  with  the 
maximum  seguency  cutoff  specified  by  the  user.  Third,  it 
calculates  the  inverse  Walsh  transform  of  the  filtered 
spectrum,  yielding  a  modified  image.  Each  of  these  steps  is 
described  in  detail  below.  The  reader  may  also  wish  to 
consult  Appendix  L,  which  contains  a  listing  of  program 
IMAGE. for. 

1.  Production  of  the  Image's  Sequeney  Spectrum 

The  program  IMAGE. for  uses  the  same  fast  Walsh  transform 
subroutine  (see  subroutine  FWT  in  Appendix  L)  as  was 
incorporated  into  program  WMASK.for,  described  earlier  in 
section  III.B.  of  this  thesis.  Two  different  types  of 
transform  geometries  were  used  in  this  research.  For  most  of 
the  images  of  interest,  a  straightforward  two-dimensional 
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Figure  7  Cutters 12 
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Figure  10  Fences 12 
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Figure  11  Emblems 12 
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resolution  degradation  was  required,  with  the  same  sequency 
filter  cutoff  in  both  X  and  Y  directions.  A  second  approach, 
which  was  applied  only  to  the  Cutter  image,  investigated  the 
effects  of  one-dimensional  sequency  filtering  along  both  the 
X  and  Y  directions,  independently,  while  leaving  the  other 
directions  unfiltered. 

The  one-dimensional  transform  procedure  is  the  simpler  of 
the  two.  If  resolution  degradation  is  desired  in  the  X 
(horizontal)  direction,  the  image  is  first  partioned  into 
rows.  The  program  calls  the  FWT  subroutine  for  that  row,  and 
returns  its  sequency  spectrum.  It  then  cycles  to  the  next 
row,  and  loops  until  all  rows'  spectra  have  been  computed. 
The  one-dimensional  transform  for  resolution  degradation  in 
the  Y  (vertical)  direction  is  effectively  the  same,  except 
that  the  program  operates  on  colvunns,  rather  than  rows. 

The  two-dimensional  transform  operation  is  essentially  a 
combination  of  both  one-dimensional  versions.  The  program 
first  performs  the  Walsh  transform  on  each  column  of  the  image 
in  turn,  then  it  calculates  the  transform  for  each  row.  This 
yields  the  two-dimensional  sequency  spectrum  of  the  original 
spectrum.  Clearly,  the  efficient  implementation  of  such  a 
numerically  intensive  calculation  is  completely  dependent  upon 
a  fast  transform  algorithm;  direct  calculation  would  be 
prohibitively  slow.  Even  with  the  FWT,  full  two-dimensional 
Walsh  spectrum  calculations  for  a  512  by  512  pixel  image 
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required  nearly  300  seconds  on  a  25  MHz  80386  PC,  using  an 
80387  math  coprocessor. 

2.  The  Sequency  Filtering  Operation 

When  the  program  IMAGE. for  is  first  invoked,  it  prompts 
the  user  for  sequency  filter  cutoff  information.  The  actual 
sequency  filtering  operation  is  then  quite  simple.  The 
program  cycles  through  each  sequency  spectrum  row  and/or 
column  that  is  to  be  filtered,  and  sets  spectral  samples  at 
sequencies  beyond  the  cutoff  sequency  to  zero.  This  is  not 
unlike  the  zeropad  procedure  described  in  section  III.B., 
although  in  this  case  the  size  of  the  spectrum  array  is  not 
increased. 

3.  Calculation  of  the  Degraded  Image 

The  inversion  of  the  filtered  sequency  spectrum  to  its 
corresponding  image  is  straightforward:  IMAGE. for  calls  FWT  to 
calculate  the  inverse  Walsh  transform  of  each  row  and/or 
column  of  the  spectrum  to  yield  the  image.  Fortunately,  this 
is  a  particularly  simple  procedure  using  the  type  of  algorithm 
developed  by  Davis  for  subroutine  FWT.  When  the  data  are 
retained  in  increasing  sequency  order  (as  they  are  here) ,  and 
the  binary  bit-reversal  sorting  portion  of  the  fast  transform 
is  followed  by  a  decimation  in  sequency,  the  Walsh  transform 
and  its  inverse  are  algorithmically  identical,  except  for  a 
simple  multiplicative  scaling  factor  following  the  inverse 


53 


transform  operation.  This  property  greatly  simplified  the 
writing  of  program  IMAGE. for 

I.  ANALYSIS  OF  RESULTS  OF  THE  RESOLUTION  DEGRADATION  PROCESS 

A  complete  set  of  resolution  degradation  results  for  each 
of  the  five  original  images  is  included  in  its  own  appendix. 
The  reason  for  this  is  that  the  conclusions  which  are  drawn 
from  these  series  of  image  photographs  are  highly  subjective; 
it  is  not  expected  that  all  readers  of  this  thesis  will  agree 
completely  with  the  interpretations.  Each  series  of  seguency 
filtered  degradation  is  represented  in  it's  entirety  for  the 
reader  to  confirm  or  refute  the  author's  analysis.  For 
illustrative  purposes,  the  image  of  the  Coast  Guard  cutter  was 
degraded  in  three  separate  manners;  once  vertically,  once 
horizontally,  and  once  in  two  dimensions.  All  other  images 
were  degraded  in  two  dimensions  only.  The  appendices  contain 
the  following  series  of  images  and  should  be  consulted 
throughout  the  following  discussion. 

•  Appendix  C  -  Cutter (512  -  0)V 

•  Appendix  D  -  Cutter (512  -  0)H 

•  Appendix  E  -  Cutter (512  -  0) 

•  Appendix  F  -  Lhouse(5l2  -  0) 

•  Appendix  G  -  Dog (512  -  0) 

•  Appendix  H  -  Fence (512  -  0) 

•  Appendix  I  -  Emblem  (5 12  0) 
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1.  One-Dimensional  Sequency  Filtering 

a.  Vertical  Filtering  of  the  Cutter  Image 

Referring  to  the  images  of  the  cutter  in  Appendix 
C,  the  effects  of  vertical  sequency  filtering  are  clearly 
illustrated.  Figure  C-1  is  the  unfiltered  image  of  the 
cutter,  accordingly  titled  Cutter512.  In  each  successive 
figure,  the  vertical  image  resolution  is  made  coarser  by  a 
factor  of  two.  For  example.  Figure  C-2  represents  Cutter256V, 
and  is  produced  from  half  as  many  sequencies  as  Cutter512. 
Figure  C-3  represents  Cutterl28V,  which  is  produced  from  half 
as  many  sequencies  as  Cutter256V,  and  so  forth.  This  pattern 
of  resolution  reduction  (in  both  one  and  two  dimensions) 
continues  throughout  this  and  all  subsequent  series  of  image 
degradations. 

The  basic  structure  of  Cutter512  in  Figure  C-1 
indicates  four  vertically  spaced  regions  corresponding  to  the 
sky,  the  hillside,  the  cutter  and  its  moorings,  and  the  water 
in  the  foreground.  With  such  an  orientation  of  dominant 
features,  one  would  expect  that  vertical  filtering  would 
preserve  these  regions  until  only  the  lowest  vertical  sequency 
components  remained.  This,  in  fact,  is  the  case.  In  Figures 
C-1  through  C-5,  corresponding  to  vertical  resolutions  of  512, 
256,  128,  64,  and  32  pixels,  this  four-region  structure 
remains  obvious.  It  begins  to  become  unidentifiable  only  at 
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vertical  resolutions  of  16  or  less,  whereas  the  cutter  itself 
loses  recognizable  form  at  a  seguency  cutoff  of  64. 

A  noticeable  feature  of  this  and  the  following  image 
degradation  series  is  the  apparent  white  "spots'*  which  seem  to 
outline  the  edges  of  the  image.  Cutter256  in  Figure  C-2 
exhibits  this  behavior  quite  clearly.  It  is  hypothesized  that 
these  spots  are  the  result  of  the  type  of  "ringing”,  or 
overshoot  phenomenon,  associated  with  the  processing  of 
adjacent  areas  of  high  contrast.  If  this  is  the  correct 
explanation  for  their  presence,  then  there  are  more  than 
likely  just  as  many  black  "spots";  the  eye  just  does  not  see 
them  standing  out  as  clearly.  The  nature  of  these  spots  and 
their  effect  on  the  imaging  process  are  offered  as  recommended 
areas  of  future  research. 

b.  Horizontal  Filtering  of  the  Cutter  Iwtage 

Appendix  D  contains  the  horizontally  degraded 
series  of  images  of  the  cutter.  As  before,  the  first  image  in 
the  series,  Cutter512  in  Figure  D-1,  represents  the  image  in 
its  undegraded  form.  All  subsequent  images  in  this  series  are 
degraded  by  a  factor  of  two  in  the  horizontal  direction. 

Cutter512  displays  much  greater  spatial  variation 
in  the  horizontal  as  compared  to  the  four  dominant  vertical 
regions  discussed  above.  Therefore,  one  would  expect  that 
eliminating  the  higher  order  sequencies  in  the  horizontal 
would  have  less  effect  on  the  overall  image  degradation.  This 
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expectation  is  validated  by  comparing  Cutter256H  in  Figiire  D-2 
with  the  original,  undegraded  image.  When  making  this  saune 
comparison  with  Cutter256V  (Figure  C-2) ,  it  becomes  readily 
apparent  that  horizontal  filtering  detracts  less  from  the 
image  quality  (as  opposed  to  vertical  filtering)  for  this 
particular  image,  due  to  it's  spatial  characteristics. 

Another  easily  seen  result  of  horizontal  filtering 
of  the  cutter  image  is  the  rapid  degradation  of  the  narrow, 
vertical  features.  The  cutter's  mast,  for  exeunple,  may  be 
described  as  such,  and  is  produced  from  higher  order 
seguencies.  One  would  then  expect  that  the  mast  would  quickly 
lose  its  detail  while  the  image  underwent  horizontal  sequency 
filtering.  This  is  indeed  the  case.  The  mast  becomes  largely 
indistinguishable  in  Cutterl28H  (Figure  D>3),  while  the 
overall  image  of  the  ship  is  still  recognizable  in  Cutter32H 
(Figure  D-5) . 

2.  Two-Dimensional  Sequenoy  Filtering 
a.  The  Coast  Guard  Cutter 

Here,  and  for  all  subsequent  image  series, 
resolution  is  diminished  by  a  factor  of  two  in  both  vertical 
and  horizontal  directions  for  each  successive  degradation.  The 
two-dimensional  filtering  of  this  image  is  displayed  in 
Appendix  E,  and  is  essentially  the  combination  of  the  one¬ 
dimensional  horizontal  filtering  as  described  above  and  the 
one-dimensional  vertical  filtering.  Analysis  of  this  series 
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concludes  that  although  the  four  distinct  vertical  regions 
(the  sky,  hillside,  moorings,  and  water)  are  still  visible 
beyond  a  sequency  cutoff  of  64,  the  spatial  structure  of  the 
cutter  itself  in  Figure  E-4  is  no  longer  recognizable. 

b.  The  Lighthouse 

The  series  of  degradations  of  this  image  appear  in 
Appendix  F.  The  undegraded  image,  Lhouse512  in  Figure  F-1, 
seems  partitioned  into  two  dominant  vertical  regions,  and 
three  horizontal.  However,  none  of  these  regions  incorporate 
any  fine  detail  in  either  dimension.  Therefore,  one  would 
expect  that  the  image  pattern  would  remain  recognizable  even 
with  comparatively  few  low-order  seguencies.  Lhouse64  in 
Figure  F-4  confirms  this,  as  the  lighthouse  is  clearly 
recognizable  when  formed  from  64  resolution  elements.  In 
fact,  Lhouse32  in  Figure  F-5,  is  still  largely  identifiable 
from  only  half  again  as  many  elements. 

c.  The  German  Shepherd 

This  series  of  degradations  appears  in  Appendix  6. 
The  undegraded  image  in  Figure  G-1  indicates  that  there  are 
few,  if  any,  right  angles  in  the  scene.  This  figure  also 
shows  a  great  amount  of  fine  detail,  as  well  as  a  somewhat 
circular  symmetry  around  the  dog's  head.  Therefore,  one  would 
expect  that  the  pattern  would  require  most  of  the  higher  order 
sequencies  to  remain  identifiable.  This  is  not  so  !  Facial 
characteristics,  such  as  eyes,  nose,  and  mouth,  appear 


58 


discernible  even  with  resolution  degraded  to  32  sec[uencies  as 
seen  in  Dog32  in  Figure  G-5.  There  is  most  likely  a  strong 
psychological  reason  for  this.  Perhaps  the  visual  information 
processing  areas  of  the  human  brain  selectively  register  this 
pattern  as  "wolf!”  or  "danger!”,  even  with  very  limited 
spatial  image  information. 

d.  The  House  with  the  Picket  Fence 

This  image  was  chosen  for  analysis  because  of  the 
pronounced  periodic  structure  exhibited  by  both  the  fence  and 
the  support  columns  on  the  front  porch.  This  is  evident  from 
Figure  H-1.  Both  features  are  vertical  structures  of  high 
contrast,  and  should  be  susceptible  to  modification  by 
horizontal  seguency  filtering.  The  spatial  period  of  the 
coltunns  is  greater  than  that  of  the  fence.  Therefore,  one 
would  predict  that  the  fence  detail  would  be  filtered  out 
first.  Fence256  and  Fencel28,  in  Figures  H-2  and  H'‘3, 
respectively,  bear  this  prediction  out.  Nonetheless,  the 
object  of  primary  concern  with  regard  to  pattern  recognition 
is  the  house  itself,  and  it  loses  recognizable  form  beyond 
resolution  with  64  elements. 

e.  The  Coast  Guard  Emblem 

This  image  was  selected  for  analysis  because  it 
represents  several  types  of  symmetry  whose  degradation  effects 
are  difficult  to  predict.  The  undegraded  image  in  Figure  I-l 
shows  some  horizontally  varying  structure  in  the  bars  of  the 
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shield,  but  the  overall  symmetry  is  circular.  Otherwise, 
there  is  very  little  purely  horizontal  or  purely  vertical 
detail,  so  that  one  is  left  to  empirical  judgment  to  evaluate 
results. 

Emblem64  in  Figure  I-*4  contributes  to  the 
conclusion  that  the  shield  structure  becomes  unrecognizable 
beyond  a  seguency  cutoff  of  64.  However,  the  bulk  circular 
outline  of  the  anchors  and  lifering  are  still  apparent;  in 
fact,  they  remain  quite  visible  in  Emblem32  in  Figure  1-5. 
This  is  somewhat  surprising  because  these  high  contrast,  high 
seguency  features  do  not  suffer  rapid  image  degradation  while 
undergoing  seguency  filtering.  As  was  the  case  with  the  image 
of  the  dog,  results  did  not  agree  closely  with  predictions. 
The  common  feature  between  these  two  departures  from 
prediction  is  the  circular  symmetry.  Recall  that  the  Walsh 
functions  which  were  used  throughout  this  project  were 
generated  in  rectangular  coordinates.  Perhaps  there  is  a 
connection  to  the  coordinate  system  of  the  applied  Walsh 
transform  and  the  dominant  symmetry  of  the  observed  scene. 
Further  research  into  this  apparent  consequence  is  in  order. 

J.  DETERMINING  THE  PATTERN  RECOGNITION  RESOLUTION  LIMIT 

From  the  seguency  filtered  degradations  of  the  five 
images,  the  author  concludes  that,  with  the  exception  of  the 
Emblem  series,  the  images  require  at  least  64  seguencies  for 
minimum  pattern  recognition.  Again,  it  must  emphasized  that 
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this  is  a  somewhat  subjective  judgment.  The  Emblem  series  is 
excluded  from  this  general  statement  because,  as  the  reader 
should  concur,  it  remains  distinguishable  beyond  this  sequency 
cutoff,  down  to  a  minimum  of  32.  It  is  believed  that  this  is 
related  to  the  highly  circular  nature  of  this  particular 
image . 

Nonetheless,  it  is  with  confidence  that  this  phase 
of  the  research  concludes  that  the  minimum  number  of 
seguencies  required  for  threshold  pattern  recognition  is  64. 
This  result  allows  one  to  quantify  several  physical 
parameters.  For  example,  this  sequency  limit  mandates  that  64 
sequencies  be  available  to  form  the  basis  set  of  orthogonal 
functions  in  each  of  the  X  and  Y  spatial  directions. 
Furthermore,  the  number  of  pixels  required  to  display  the 
output  image  must  be  fixed  at  4096,  which  represents  the 
square  of  64. 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 


The  intent  of  the  initial  task  of  this  project  was  to 
mathematically  model  the  two-dimensional  sequency-ordered 
Walsh  functions,  and  to  determine  the  diffraction  pattern  of 
each.  This  was  accomplished.  The  result  of  studying  the 
patterns  lead  to  the  formulation  of  the  following  equation 
which  quantifies  the  diffraction  limited  performance  of  the 
encoding  scheme: 

1.^  sindi  <  6  . 

The  angle  6  represents  the  maximum  tolerable  diffraction  angle 
which  will  ensure  that  all  of  the  diffracted  photon  flux 
passing  through  the  Walsh  mask  will  be  focused  onto  the 
detector.  This  angle,  coupled  with  the  size  of  the  detector 
and  the  wavelength  of  operation,  provides  useful  information 
to  the  follow-on  work  which  will  proceed  to  design  the  optics 
for  a  prototype  instrument.  This  simple  relationship  is 
illustrated  in  Figure  6. 

The  second  question  addressed  by  this  thesis  was  directed 
towards  establishing  the  minimum  number  of  sequencies  required 
to  produce  a  recognizable  image.  This  was  accomplished  by 
generating  images  with  increasing  numbers  of  higher  order 
sequencies  filtered  out.  Once  the  image  pattern  could  no 
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longer  be  deemed  identifiable,  it  was  concluded  that  the 
minimum  pattern  recognition  limit  had  been  reached.  From 
examining  the  all  of  the  series  of  seguency  filtered  images  in 
Appendices  E  -  I,  it  was  concluded  that  no  less  than  64 
seguencies  are  reguired  to  produce  a  recognizable  pattern. 

It  is  recommended  that  future  thesis  work  be  directed 
toward  investigating  how  the  coordinate  system  of  the  Walsh 
functions  affect  the  results.  This  guestion  is  raised  from 
analysis  of  the  series  of  degraded  images  of  the  circular 
Coast  Guard  emblem  (Appendix  I) .  It  appears  that  perhaps  the 
dominant  symmetry  of  the  object  plays  an  important  role  in 
determining  the  minimum  number  of  seguencies  reguired  to 
resolve  it.  An  interesting  project  may  be  to  regenerate  the 
findings  of  this  project  while  using  Walsh  functions  derived 
under  a  different  coordinate  system,  and  make  comparisons 
between  the  two. 

Furthermore,  gueries  into  the  dependency  of  image 
resolution  as  a  function  of  the  available  greyscale  level,  and 
the  cause  and  effects  of  the  "ringing"  phenomenon  apparently 
caused  by  processing  adjacent  elements  of  high  contrast,  are 
recommended.  Either  aspect  could  possibly  affect  the  findings 
of  this  thesis,  and  may  mislead  future  research  based  upon  the 
conclusions  drawn  herein. 
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APPEMDIZ  A 


«  LT  B  .J.  nUSSELMAN 

« 

*  THESIS  PROGRAM  -  MMASK.FOR 

* 

*  THIS  PROGRAM  COMPUTES  MATHEMATICAL  MODELS  OF  WALSH  SEOUENCY 

*  MASKS,  AND  CALCULATES  THEIR  DIFFRACTION  PATTERNS  THROUGH 

*  2  -  D  FOURIER  TRANSFORMS.  THE  COMPLEX  TRANSFORM  MATRIX  IS 

*  REDUCED  TO  AN  INTENSITY  MATRIX,  WHICH  IS  STORED  IN  CARTESIAN 

*  COORDINATES  <X,Y,Z)  IN  A  FILE  CALLED  Ct\FORTRAN\TEMP\XYZ.DAT 

* 

«  THE  INPUT  MATRIX  CWALSH  MASK)  IS  16  X  16,  AND  IS  CENTERED 

*  IN  A  32  X  32  ZEROPAD.  THE  ENTIRE  MATRIX  IS  SENT  TO  THE 

*  2  >  D  FOURIER  TRANSFORM  SUBROUTINE. 

* 

INTEGER  ISIGN,NDIM, ROW , NROW , COL , NCOL , NN , N , R , J , SEQ , NUMS AMPLS 
DIMENSION  NN  (2) 

REAL  I DATA (64,32), TDATA < 2048) , PDATA ( 1 024 ) 

REAL  W(64) , X (961 ) , Y(961 ) , Z (961 ) , TEMP,EESTAR(32, 32) 

REAL  A(e,8),B(8, 16),C(16, 16) ,CLIN(256) 

NDIM  «  2 

• 

•  NN  (1)  REPRESENTS  ROWS,  NN  (2)  REPRESENTS  COLUMNS 

* 

NN(1)  -  32 
NN(2)  >  32 

« 

NCOL  -  32 
NROW  s  32 

* 

ISIGN  -  1 
NUMSAMPLS  -  64 

* 

OPEN  (11, FILE  -  ’Ci\FORTRAN\TEMP\XYZ.DAT’,  STATUS  *  »OLD») 

• 

WRITE(*,*)  'WMASK.FOR  IS  CURRENTLY  CONFIGURED  TO  PRODUCE  A  BASIS' 
WRITE(*,*)  'SET  OF  64  WALSH  SEQUENCIES  NUMBERED  FROM  O  TO  63.' 
WRITE(*,*) 

WRITE(*,*)  'ENTER  SEOUENCY  NUMBER  TO  GENERATE  AND  TRANSFORM  :' 
REAO(*,*)  SEQ 

* 

DO  S2  ROW  >  1,  NROW  •  2 
DO  32  COL  -  1,  NCOL 

I DATA (ROW, COL)  -  0.0 
52  CONTINUE 

* 

» 

00  10  J  «  1,  NUMSAMPLS 
10  W(J)  -  0.0 
W(SEQ  *  1)  -  1.0 

• 

CALL  FWT  (NUMSAMPLS, W, .FALSE. ) 

• 

♦  THIS  SECTION  DEVELOPS  THE  2X  EXPANDED  MATRIX  USING  WALSH 

*  SEQUENCIES  AS  RETURNED  FROM  SUBROUTINE  FWT.  THE  MATRIX 

*  IS  WRITTEN  TO  THE  SCREEN  AS  A  VISUAL  CHECK. 

• 

J  ■  1 

DO  32  ROW  «  1,  SORT (NUMSAMPLS) 

DO  32  COL  «  1,  SORT (NUMSAMPLS) 
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A(ROW,COL)  -  .5  *  (W(J>  ♦  1) 

J  -  J  +  1 

32  CONTINUE 

DO  43  ROW  -  1,  SaRT(NUMSAHPLS) 

DO  43  COL  -  1,  SQRT(NUMSAnPLS> 

B<ROW,  2  «  COL>  «  A(ROW,COL> 

B(ROW,  2  *  COL  -  1>  >  B(ROU,  2  «  COL) 

43  CONTINUE 

DO  34  ROM  -  1,  SQRT(NUnSAMPLS) 

DO  34  COL  >  1,  2  «  SQRTINUnSAMPLS) 

C(2  «  ROW  ,  COL)  «  B(RQW,COL) 

C(2  *  ROU  -  l.COL)  -  C<2  *  ROW, COL) 

34  CONTINUE 

DO  3S  ROW  «  1,  2  «  SQRT(NUnSAMPLS) 

WRITE<*,201)  <C<ROW,COL)  ,  COL  -1,2*  SQRT(NUMSAMPLS) ) 

35  CONTINUE 

READ(*,*) 

J  -  1 

DO  36  ROW  -  1,  2  «  SQRT(NUnSAnPLS) 

DO  36  COL  «  1,  2  •  SORTCNUliSAMPLS) 

CLIN(J)  -  C(ROW,COL) 

J  -  J  ♦  1 

36  CONTINUE 


*  TORN  THE  COMPLEX  MATRIX  WITH  ALL  IMAGINARY  COMPONENTS  ZERO 

« 

CALL  ZEROPAD  (CLIN, IDATA) 

« 

CALL  rLIPFLOP(IDATA,NROW,NCQL) 


CALL  MAT2L1N  <NROW,NCOL, IDATA, TDATA) 
CALL  FOURN  (TDATA, NN, NDIM, ISIGN) 

« 

CALL  L1N2MAT  (NROW, NCOL, IDATA, TDATA) 

« 

CALL  FLIPFLOP  ( IDATA, NROW, NCOL) 

* 

CALL  INTENSITY  (NROW, NCOL, IDATA, EESTAR) 

* 

CALL  MAT2XYZ  (NROW, NCOL, EESTAR, X,Y,Z) 

» 

*  WRITE  XYZ  MATRIX  TO  FILE  FOR  PLOTTING 

* 

DO  15  N  -  1,  (NROW  -  1)  ♦  (NCOL  -  1) 
WRITEdl,  102)  X(N),  Y(N),  Z(N) 

15  CONTINUE 

* 

100  FORMAT (1X,F5. 2) 

101  FORMAT ( IX, 8 (Fa. 2, 2X)) 

102  F0RMAT(1X,F8.3,2X,F8.3,2X,E12.6) 

201  FORMATdX,  16(F4.  1,  IX)) 

END 


» 

SUBROUTINE  ZEROPAD  (CLIN, IDATA) 


INTEGER  J, ROW, COL 

REAL  CLINd024)  ,  IDATA(P4,32) 
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J  ■  1 

DO  12  ROW  -  17,47,2 
DO  12  COL  ■  9,24 

I DATA (ROW, COL)  «  CLINCJ) 

J  -  J  1 
12  CONTINUE 
RETURN 
END 

SUBROUTINE  MAT2LIN  (NROW, NCOL, IDATA, TDATA) 

INTEGER  N, ROW, COL 

REAL  I DATA (64, 32), TDATA (2048) 

N  *  1 

DO  4  COL  *  1,  NCOL 

DO  4  ROW  «  1,  NROW  *  2 

TOATA(N)  -  I DAT A (ROW, COL) 

N  -  N  +1 
4  CONTINUE 
RETURN 
END 

SUBROUTINE  LIN2f1AT  (NROW, NCOL,  IDATA,  TDATA) 

INTEGER  N, COL, ROW 

REAL  I DATA (64, 32), TDATA (2048) 

N  *  1 

DO  6  COL  -  1,  NCOL 

DO  6  ROW  «  1,  NROW  *  2 
I OAT A (ROW, COL)  ■  TOATA(N) 

N  •  N  +  1 
6  CONTINUE 
RETURN 
END 

SUBROUTINE  INTENSITY  (NROW, NCOL, IDATA, EESTAR) 

INTEGER  N, ROW, COL 

REAL  I0ATA(64, 32) ,PDATA( 1024) , EESTAR (32, 32) 

N  ■  1 

DO  8  COL  ■  1,  NCOL 

DO  8  ROW  «  1,  (2  «  NROW)  -1,  2 

PDATA(N)  -  (IDATA(R0W,C0L))»*2  ♦  (IDATA(ROW  +  1,C0L))**2 
N  «  N  ♦  1 

8  CONTINUE 

N  -  1 

DO  9  COL  -  1,  NCOL 
DO  9  ROW  «  1,  NROW 

EESTAR (ROW, COL)  -  PDATA(N) 

N  ■  N  1 

9  CONTINUE 
RETURN 
END 

SUBROUTINE  nAT2XYZ  (NROW, NCOL, EESTAR, X, Y, Z ) 
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integer  N, ROM, col 

REAL  EESTARC32,32) ,X<961>,Y<S61>,Z<961> 


N  »  1 

DO  2  ROW  «  1,  NROU  -  1 
DO  2  COL  •  1,  NCOL  -  1 

X<N)  «  -6.2831  <ROW  -  1>  *  (12.5&62/30> 
Y<N>  -  -6.2831  *  CCOL  -  1)  ♦  <12. 5662/30) 
Z<N)  -  EESTARCROW  ♦  1,C0L  ♦  1) 

N  ■  N  ♦  1 
CONTINUE 
RETURN 
END 

SUBROUTINE  FLIPFLOR  (BOX, NROW, NCOL) 

INTEGER  R,C, NROW, NCOL 
REAL  BOX (64, 32) 

RE-ORDER  MATRIX  LEFT  TO  RIGHT 

DO  2  C  -  1,  NCOL/2 

DO  2  R  -  1,  NROW  *  2 
TEMP  •  BOX(R,C) 

BOX(R,C)  -  BOX  <R,  NCOL/2  C) 

BOX  (R,  NCOL/2  C)  •  TEMP 

CONTINUE 

RE  -  ORDER  MATRIX  TOP  TO  BOTTOM 

DO  3  C  •  1,  NCOL 

DO  3  R  •  1,  <NR0W*2)/2 
TEMP  -  BOX<B,C) 

BOX<R,C)  -  B0X((NROW*2)/2  ♦  R,C) 
B0X((NR0W*2)/2  ♦  R,C)  -  TEMP 
CONTINUE 

return 

END 

SUBROUTINE  FOURN  (DATA,NN,NDIM, ISIGN) 

REAL  UR , W I , WPR , WP I , WTEMP , THETA 
DIMENSION  NN(NDIM),DATA(*) 

NTOT  -  1 

DO  11  IDIM  «  1,  NDIM 

NTOT  ■  NTOT  *  NNdDIM) 

CONTINUE 
NPREV  «  1 

DO  18  IDIM  -  1,  NDIM 
N  -  NNdDIM) 

NREM  »  NTOT  /  (N  *  NPREV) 

IPl  •  2  *  NPREV 
IP2  *  IPl  *  N 
1P3  -  1P2  *  NREM 
12REV  •  1 

DO  14  12  ■  l,IP2,IPt 

IF  (12  .LT.  I2REV)  THEN 

DO  13  II  •  12,  12  ♦  IPl  -  2,  2 


o  o  n  n  o  o 


DO  12  13  >  II,  1P3,  1P2 

I3REV  -  I2REV  +13-12 
TEMPR  -  DATA  (13) 

TEMPI  -  DATA  (13  ♦  1) 

DATA  <I3>  ■  DATA  (I3REV) 

DATA  (13  +  1)  •  DATA  (I3REV  ♦  1) 

DATA  (13REV)  -  TEMPR 
DATA  (13REV  ♦  1 )  •=  TEMPI 

12  CONTINUE 

13  CONTINUE 
END  IF 

IBIT  -  IP2/2 

1  IF  ((IBIT  .QE.  IPl)  .AND.  (I2REV  .GT.  IBIT))  THEN 

I2REV  »  I2REV  -  IBIT 
IBIT  «  IBlT/2 
GO  TO  1 
END  IF 

I2REV  -  I2REV  ♦  IBIT 

14  CONTINUE 
IFPl  -  IPl 

2  IF  (IFPl  .LT.  1P2)  THEN 

1FP2  ■  2  *  IFPl 

THETA  »  ISIQN  *  6.28318S30717179S9D0/  (1FP2/IP1) 
gpR  -  -2. DO  *  DSIN(O.SDO  *  THETA)  **  2 
WPI  -  DSIN  (THETA) 

MR  -  l.DO 
MI  >  O.DO 

DO  17  13  -1,  IFPl,  IPl 

DO  16  II  -  13,  13  +1P1  -  2,  2 
DO  15  12  -  11,  XP3,  1FP2 
K1  -  12 
K2  ■  K1  ♦  IFPl 

TEMPR  >  8NGL(UR)*DATA(K2)-SNGL(WI)*DATA(K2+1) 
TEMPI  »  SNSL(MR)*DATA(K2+l>+SNGl.(Ml)*DATA(K2) 
0ATA(K2)  -  DATA(Kl)  -  TEMPR 
0ATA(K2  ♦  1)  ■  DATAdCl  +  1)  -  TEMPI 
OATA(Kl)  •  DATA(Kl)  ♦  TEMPR 
DATA(Kt  ♦  1)  ■  DATA(K1  ♦  1)  ♦  TEMPI 

15  CONTINUE 

16  CONTINUE 
WTEMP  ■  MR 

UR  -  UR  *  UPR  -  U1  *  UPI  ♦  UR 
UI  -  UI  ♦  UPR  ♦  UTEMP  *  UPI  +  UI 

17  CONTINUE 
IFPl  «  IFP2 
GO  TO  2 
END  IF 

NPREV  »  N  *  NPREV 

18  CONTINUE 
RETURN 
END 

SUBROUTINE  FUT  (NUM8AMPLS, DAT ARRAY, ENABLESCALE) 

FAST  UALSH  TRANSFORM  (TRANSLATED  FROM  D.S.D.’S  VERSION  IN  “C”) 

NUMSAMPLS  e  NUMBER  OF  SEQUENCIES  AND  SAMPLING  INTERVALS 
IN  UALSH  BASIS  (INTEGER) 

DATARRAY  *  INPUT/OUTPUT  ARRAY  (REAL) 
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o  o  o 


ENABLESCALE  -  TRANSFORM  SCALING  ENABLE  SWITCH 

INTEGER  NUMSAMPLS 
REAL  DATARRAYCn 
LOGICAL  ENABLESCALE 

CALL  BITREVSORT  (NUMSAMPLS, DAT ARRAY) 

CALL  DECINSEQ  (NUMSAMPLS, DAT ARRAY) 

IF  (ENABLESCALE)  CALL  SCALE  (NUMSAMPLS, DATARRAY) 

RETURN 
END 

SUBROUTINE  BITREVSORT  (NUMSAMPLS, DATARRAY) 

C  BINARY  BIT  -  REVERSAL  SORT 

» 

I NTEGER  NUMSAMPLS , NUMSD2 , D I RECT , REVERSED , DAT AD , DAT AR , OFFSET 
REAL  DATARRAY ( 1 ) , SCRATCH 

« 

NUMSD2  *  NUMSAMPLS  /  2 
REVERSED  -  O 

DO  30  DIRECT  ■  O,  NUMSAMPLS  -  1 

IF  (DIRECT  .GE.  REVERSED)  GO  TO  10 
OATAD  mi*  DIRECT 
DATAR  -  1  REVERSED 
SCRATCH  *  DATARRAY (DATAD) 

DATARRAY (DATAD)  *  DATARRAY (DATAR) 

DATARRAY (DATAR)  >  SCRATCH 
10  OFFSET  ■  NUMSD2 

20  IF ((OFFSET  .6T.  REVERSED)  .OR.  (OFFSET  .LT.  2))  GO  TO  30 
REVERSED  -  REVERSED  -  OFFSET 
OFFSET  -  OFFSET  /  2 
GO  TO  20 

30  REVERSED  >  REVERSED  *  OFFSET 
RETURN 
END 

* 

SUBROUTINE  DECINSEQ  (NUMSAMPLS, DATARRAY) 

* 

C  DECIMATION  IN  SEQUENCE 

* 

INTEGER  NUSAMPLS, L0G2NUMSAMPLS, MSB, BLOCKS I ZE, STAGE 
I NTEGER  NODE 1 START , NaOE2START , NODE 1 , N0DE2 
LOGICAL  TOGGLE 
REAL  DATARRAY ( 1 ) , SCRATCH 

« 

L0G2NUMSAMPLS  «  O 
MSB  -  NUMSAMPLS 
10  MSB  «  MSB  /  2 

IF  (MSB  .EQ.  0)  GO  TO  20 

L0G2NUMSAMPLS  -  L0G2NUMSAMPLS  *  1 
GO  TO  10 

20  BLOCKS IZE  >  NUMSAMPLS  /  2 

DO  80  STAGE  >  1,  L0G2NUMSAMPLS 
TOGGLE  -  .FALSE. 

NODE 1 START  -  0 

30  IF  ( NODE 1 START  .GE.  NUMSAMPLS)  GO  TO  BO 


69 


* 


* 

c 

* 


« 


NQDE2START  -  NODEISTART  ♦  BLOCKSI2E 
N0DE2  «  N0DE2START 
IF  (.NOT.  TOOQLE)  80  TO  50 

DO  40  NOOEl  ■  NODEISTART, N0DE2ST ART  -  1 
SCRATCH  ■  DATARRAY(N0DE1  1) 

DAT  ARRAY  C  NODE  t*l) -DAT ARRAY ( NODE 1 + 1 ) -DATARRA Y  <  N0DE2+ 1 ) 
DATARRAY(N0DE2-^1>  -  DATARRAY(N0DE2-^1  >  ♦  SCRATCH 
N0DE2  *  N0DE2  *  1 
GO  TO  70 

DO  60  NODEl  -  NODEISTART, N0DE2ST ART  -  1 
SCRATCH  -  DATARRAY(NODEl-*-l) 

DATARRAY(N0DE1+1)-DATARRAY(N0DE1+1 )+DATARRAY(N0DE2+l ) 

DATARRAY(N0DE2-«-l>  -  SCRATCH  -  DATARRAY(N0DE24-1  ) 

N0DE2  -  N0DE2  *  1 
TOGGLE  -  .NOT.  TOGGLE 
NODEISTART  -  N0DE2 
GO  TO  30 

BO  BLOCKS IZE  »  BLOCKS I ZE  /  2 
RETURN 
END 


40 


50 


60 

70 


SUBROUTINE  SCALE  (NUMSANPLS, DAT ARRAY) 

UALSH  TRANSFORM  SCALING 

INTEGER  NUMSAMPLS, I 

REAL  DA TARRAY(1),SCALEF ACTOR 

SCALEFACTOR  -  1.0  /NUMSAMPLS 
DO  10  I  -  1,  NUMSAMPLS 

DATARRAY(l)  •  DATARRAYCI)  *  SCALEFACTOR 
10  CONTINUE 
RETURN 
CLOSE  (11) 

END 
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n^ure  B-«  Diffraction  Pattern  for  Sequency  55 
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APPENDIX  C 


Figure  C-l  Cutter5l2V 
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Figure  c-2  Cutter256V 
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Figure  C-3  Cutter 128V 
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Figure  C-4  Cutter64V 
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Figure  C-5  Cutters 2V 
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Figure  C-6  Cutterie 
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Figiire  C-8  Cutter4V 
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APPENDIX  D 


Figure  D-1  Cutter512H 
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Figure  D-2  Cutter 2 56H 
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Figure  D-3  Cutter 12 8H 
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Figure  D-4  Cutter64H 


Figure  D-5  Cutter 3 2H 


Figure  D-6  Cutter 16H 
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Figure  D-7  CutterSH 
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Figure  D-8  Cutter4H 
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Figure  D-9  Cutter2H 


Figure  D-10  CutterlH 
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Figure  E~1  Cutters 12 
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Figure  E-2  Cutter 2 56 
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Figure  E-4  Cutter64 
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Figxire  E-5  Cutter 3  2 
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Figure  E-6  Cutter 16 
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Figure  E-10  Cutterl 
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APPENDIX  F 


Figure  F-1  Lhouse512 
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Figure  F-2  lihou8e256 
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Figure  F-3  Lhousel28 
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Figure  G-l  Dog512 
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Figure  G-3  D09I28 
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Figure  H-1  Fence512 
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Figure  H-2  Fence256 
148 


Figiire  H-4  Fence64 
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APPEMDIZ  I 


Figure  I-1  Einblein512 
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Figure  1-2  Eittblem256 
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APPENDIX  J 


/********«**************•**< 

/*  C512:  Conversion  of  512  X  512  OT-IRIS  Inages  into  Sinple  Byte  Fomat  •/ 

/"  910614DSD  Microsoft  C5.1  */ 

y ************************************************************************** ! 


/include  <stdio.h> 
/include  <stdlib.h> 


nain  (void) 

{  unsigned  int  a,  row,  column; 

char  instreanname [40],  outstreannane [40]; 

FILE  *  instream,  *  outstream; 

system  ("CLS"); 

for  (a  *  1;  a  <*  79;  a  ++) 

Dutch  ('*'); 
printf  ("\n\r"); 
for  (a  »  1;  a  <*  17;  a  ++) 

Dutch  ( '  ' ) ; 

printf  ("T512:  512  X  512  DT-IRIS  Inage  Conv6rsion\n\r") ; 
for  (a  *  1;  a  <»  79;  a  ++) 
putch  ('*'); 
printf  ("\n\n\r"); 

printf  ("Name  of  DT-IRIS  input  file  ?  "); 
scant  ("ts",  instreamnane) ; 

if  ((instream  *  fopen  (instreanname,  "rb"))  NULL) 

{  fprintf  (stderr,  "Can't  open  input  file!\n"); 
exit  (0) ;  } 
printf  ("\n\n\r") ; 

printf  ("Name  of  translated  output  file  ?  "); 
scanf  ("ts",  outstreannane); 

if  ((outstream  ■  fopen  (outstreannane,  "wb"))  *■  NULL) 

{  fprintf  (stderr,  "Can't  open  output  file!\n"); 
exit  (0) ;  } 
printf  ("\n\n\r"); 
for  (a  ■  1;  a  <»  512;  a  ++) 
f getc  ( instream) ; 
for  (row  ■  1;  row  <■  512;  row  ++) 

{  for  (column  ■  1;  column  <■  512;  column  ++) 

{  a  f getc  ( instream) ; 

fputc  (a,  outstream);  } 
fputc  ('\n',  outstream);  } 
f close  (instream) ; 
f close  (outstream) ; 

printf  ("Output  file  contains  512  lines  of  512  single*byte  unsigned  ") ; 
printf  ("integers\n\n\n\n\r") ; 
return;  > 


/ 
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APPENDIX  K 


^ ***************************************************************** ******* 

/*  D512:  VGA  Display  for  512  X  512  Single-byte  Greyscale  Images 

/*  910614DSD  Microsoft  C5.1 

/************************************************************************* 

#include  <stdio.h> 

^include  <stdlib.h> 

/include  <conio.h> 

/include  <graph.h> 

/include  <math.h> 

main  (void) 

(  int  basis,  a,  amin,  amax,  row,  column,  rowoffset,  xstart,  y,  x,  color; 
long  int  tints [16]; 
char  instreamname [40]; 

FILE  *  instream; 
struct  videoconfig  vc; 
basis  >  512; 
tints[0]  »  0x00000000; 
tints[l]  ■  0x00040404; 
tints [2]  -  0X00080808; 
tints[3]  -  OxOOOcOcOc; 
tints[4]  «  0x00101010; 
tints [5]  -  0x00141414; 
tints [6]  -  0x00181818; 
tints[7]  ■  OxOOlclclc; 
tints[8]  -  0X00202020; 
tints[9]  >  0x00242424; 
tints[10]  «  0x00282828; 
tints [11]  ■  Cx002c2c2c; 
tints[12]  ■  0x00303030; 
tints[13]  -  0x00343434; 
tints[14]  >  0x00383838; 
tint>:.[15]  ■  0x003c3c3c; 
system  ("CLS"); 
for  (a  ■  1;  a  <■  79;  a  ++) 
putch  ('*'); 
printf  ("\n\r"); 
for  (a  «  1;  a  <■  20;  a  ++) 
putch  ( '  ' ) ; 

printf  ("0512:  512  X  512  Greyscale  Image  Display\n\r") ; 
for  (a  -  1;  a  <«  79;  a  ++) 
putch  ('*'); 
printf  ("\n\n\r"); 
printf  ("Name  of  input  file  ?  "); 
scanf  ("%s",  instreamname); 

if  ((instream  ■  fopen  (instreamname,  "rb"))  NULL) 

{  fprintf  (stderr,  "Can't  open  input  data  file!\n"); 
exit  (0);  } 
printf  ("\n\n\r"); 
amin  amax  >0; 

for  (row  “1;  row  <»  basis;  row  ++) 

{  for  (column  *  1;  column  <«  basis;  column  ++) 

(  a  K  f getc  ( instream) ; 
if  (a  <  amin) 
amin  *  a; 

else  if  (a  >  amax) 
amax  -  a;  } 
fgetc  (instream) ;  } 
rewind  (instream); 

_sutvideomode  (_VRES16COLOR) ; 

_remapallpalette  (tints) ; 

_getvideoconf ig  (8  vc) ; 

_clearscreen  (_GCLEARSCREEN) ; 
rowoffset  -  (basis  -  vc.ntunypixels)  /  2; 
for  (row  -  i;  row  <»  rowoffset;  row  ++) 
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^®tart  -  (vc.nimxpixels  -  basis)  /  2; 

for  <“  vc.numypixels;  row  ++) 

for  (column  «  l;  column  <-  basis;  column  ++) 

{  a  «  fgetc  (instream) ;  ' 

color  -  floor^  Uvc.nrmcolors  -  l)  *  (a  -  amin)  /  (amax 

_8etcolor  (color) • 

_8etpixel  (X,  y); 

++  x;  } 

fgetc  (instream); 

,  ,  y?  > 

f close  (instream); 
while  (!  Jcbhit  ()) 

_clearscreen  (  GCLEARSCREEN) ; 

_setvAdeomode  7  DEFAULTMODE) ; 
return;  }  ~ 


/ 


5 


amin)  + 
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o  n  o  o 


APPENDIX  L 


*  LT  B.J.  MUSSELMAN 

« 

C  THESIS  PROGRAM  -  IMAGE. TOR 

THIS  PROGRAM  IS  CONFIGURED  TO  PERFORM  A  2-DIMENSIONAL  WALSH 
TRANSFORM  ON  AN  INPUT  IMAGE  FILE.  THE  INVERSE  TRANSFORM  IS 
THEN  TAKEN  AND  THE  TRANSFORMED  DATA  IS  WRITTEN  TO  SEPARATE 
OUTPUT  IMAGE  FILE. 

I NTEGER*2  NUMROWS , NUMCOLS , MAXSEQ , CUTOFF , SEQUENC  Y , F I LTER ( 5 1 3 ) 
INTEGER*2  ROW, COL 
INTEGER*!  BDAT <512, 51 3) 

1NTEGER*4  W<513) 

NUMROWS  =  512 
NUMCOLS  =  513 
MAXSEQ  =  NUMROWS  -  1 

OPENCll,  FILE  *  »Bs\EMBLEM*,  FORM=»  BINARY',  STATUS  ■  'OLD') 

OPEN  <12,FILE*'C:\ STUDENTSV  BR I ANN  I MAGE ' , FORM* ' B I NARY ' , STATUS- ' OLD ' ) 

WRITE <♦,*)  'ELIMINATE  SEQUENC Y  COEFFICIENTS  BEYOND  -  ' 

READC*,*)  CUTOFF 

C  SEQUENCY  FILTER  DEFINED  IN  THIS  SECTION 

• 

DO  2  SEQ  *  0,  MAXSEQ 
FILTER<SEQ  +  1)  ■  1 

IFCSEQ  .6T.  CUTOFF)  FILTER<SEQ  ♦  1 )  -  0 

2  CONTINUE 

* 

C  READ  IMAGE  FILE,  ROW-BY-ROW  INTO  ARRAY  BDAT 

« 

DO  4  ROW  -  1,  NUMROWS 
DO  3  COL  -  1,  NUMCOLS 

READ(ll)  BDAT (ROW, COL) 

3  CONTINUE 

4  CONTINUE 

• 

C  THIS  SECTION  PERFORMS  VERTICAL  SEQUENCY  FILTERING 

• 

DO  8  COL  -  1,  NUMCOLS  -  1 
WRITE(*,*)  COL 
DO  5  ROW  =  1,  NUMROWS 

WCROWi  *  BDAT (ROW, COL 1 

5  CONTINUE 

« 

CALL  FWT  vNUMROWS,  W,  .FALSE.) 

* 

DO  6  SEQ  =  O,  MAXSEQ 

WfSEQ  ♦  l  i  -  WCSEQ  *  1)  *  FILTER(SEQ  1) 

6  CONTINUE 

• 

CALL  FWT  (NUMROWS,  W,  .TRUE.) 

« 

DO  7  ROW  =  1 ,  NUMROWS 

BDAT(ROW,COL)  -  W(ROW) 

7  continue 

e  continue 
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o**  **♦  ooooooooo*  *  *n*  ♦*  **  *o 


THIS  SECTION  PERFORMS  HORIZONTAL  SEQUENCY  FILTERING 

DO  IB  ROW  >  1,  NUMROWS 
WRITE!*,*)  ROW 
DO  15  COL  «  1,  NUMCOLS 
W(COL)  >  BDAT!RQW,COL> 

15  continue 

CALL  FWT  ((NUMCOLS  -  1),  W,  .FALSE.) 

DO  IS  SEQ  «  0,  MAXSEQ 

w(sea  *  1)  «  w(SEQ  ♦  i)  *  filter(seo  ♦  i) 

IS  CONTINUE 

CALL  FWT  ((NUMCOLS  -  1),  U,  .TRUE.) 

DO  17  COL  »  1,  NUMCOLS  -  1 
BDAT(ROW,COL>  «  W(COL) 

17  CONTINUE 

18  CONTINUE 

WRITE  IMAGE  FILE,  ROW-BY-ROW 

DO  11  ROW  -  1,  NUMROWS 
DO  13  COL  *  1,  NUMCOLS 

WRITE! 12)  BOAT (ROW, COL) 

13  CONTINUE 
11  CONTINUE 

CLOSE  (11) 

CLOSE  (12) 

STOP 

END 

SUBROUT I NE  FWT  <  NUMSAMPLS , DATARRAY , ENABLESC ALE  > 

FAST  WALSH  TRANSFORM  (TRANSLATED  FROM  D.S.D.’S  VERSION  IN  "C“> 

NUMSAMPLS  ■  NUMBER  OF  SEQUENCIES  AND  SAMPLING  INTERVALS 
IN  WALSH  BASIS  (INTEGER) 

DATARRAY  ■  INPUT/OUTPUT  ARRAY  (REAL) 

ENABLESCALE  «  TRANSFORM  SCALING  ENABLE  SWITCH 

INTE6ER*2  NUMSAMPLS 
INTEQER*4  DATARRAY(l) 

LOGICAL  ENABLESCALE 

CALL  BITREVSORT  (NUMSAMPLS, DATARRAY) 

CALL  DECINSEO  (NUMSAMPLS, DATARRAY) 

IF  (ENABLESCALE)  CALL  SCALE  (NUMSAMPLS, DATARRAY > 

RETURN 

END 

SUBROUTINE  BITREVSORT  (NUMSAMPLS, DATARRAY) 

BINARY  BIT  -  REVERSAL  SORT 
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INTEGER  NUMSD2 , D 1 RECT , REVERSED , DATAD , DATAR , OFESET 
INTEQER*2  NUMSAMPLS 
REAL  SCRATCH 
INTEGER*4  DATARRAYd) 

* 

NUMSD2  -  NUMSAMPLS  /  2 
REVERSED  «  0 

DO  30  DIRECT  •  0,  NUMSAMPLS  -  1 

IF  <DIRECT  .QE.  REVERSED)  GO  TO  10 
DAT AD  •  I  *  DIRECT 
DATAR  «  1  REVERSED 
SCRATCH  -  DAT ARRAY (DAT AD) 

DATARRAY(DATAD)  «  DATARRAYC DATAR) 

DATARRAY( DATAR)  -  SCRATCH 
10  OFFSET  «  NUMSD2 

20  1F( (OFFSET  .6T.  REVERSED)  .OR.  (OFFSET  .LT.  2))  60  TO  30 

REVERSED  ■  REVERSED  -  OFFSET 
OFFSET  ■  OFFSET  /  2 
GO  TO  20 

30  REVERSED  -  REVERSED  OFFSET 
RETURN 
END 

« 

SUBROUTINE  DECINSEQ  (NUMSAMPLS, OAT ARRAY) 

C  DECIMATION  IN  SEQUENCE 

« 

I NTEGER  L0G2NUMSAMPLS , MSB , BLOCKS 1 ZE , STAGE 
I NTEGER  NODE 1 START , N0DE2START , NODE 1 , N0DE2 
INTEGER«2  NUMSAMPLS 
INTEGER«4  DATARRAYd) 

LOGICAL  TOGGLE 
REAL  SCRATCH 

* 

L0G2NUMSAMPLS  «  0 
MSB  «  NUMSAMPLS 
10  MSB  >  MSB  /  2 

IF  (MSB  .EQ.  0)  GO  TO  20 

L062NUMSAMPLS  -  L0G2NUMSAMPLS  *  1 
GO  TO  10 

20  BLOCKS I ZE  -  NUMSAMPLS  /  2 

DO  80  STAGE  >  1,  L0G2NUMSAMPLS 
TOGGLE  -  .FALSE. 

NODE 1 start  •  0 

30  IF  ( NODE 1 start  .GE.  NUMSAMPLS)  60  TO  80 
N0DE2START  -  NODEISTART  ♦  BLOCKSIZE 
N00E2  «  N0DE2START 
IF  (.NOT.  TOGGLE)  GO  TO  50 

DO  40  NODEl  »  NODEISTART, N0DE2ST ART  -  1 
SCRATCH  *  DATARRAY (NODEl  1) 

DATARRA Y (NODE I ♦ X ) «DAT ARRAY ( NODE 1 ♦ I > -DATARRAY ( N0DE2^ 1 > 
DATARRAY  (N0DE2-*-l)  •  DATARRAY  ( N0DE2  *•  1 )  +  SCRATCH 
40  N0DE2  ■  N0DE2  1 

GO  TO  70 

50  DO  60  NODEl  •  NODEISTART, NDDE2START  -  1 

SCRATCH  ■  DATARRAY (NODEI+ I ) 

DATARRAY  (  NODE  1  ♦  1 )  -DATARRAY  (  NODE  1  X  )  ♦DATARRAY  (  N0DE2+ 1  ) 
DATARRAY (N0DE2-*-l)  «  SCRATCH  -  DATARRAY (N0DE2+1 ) 

60  N0DE2  «  N0DE2  ♦  1 

70  toggle  «  .NOT.  TOGGLE 
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NODE  1 START  *  N0DE2 
QO  TO  30 

80  BLOCK-SIZE  »  BLOCKS! ZE  /  2 
RETURN 
END 

* 

SUBROUTINE  SCALE  (NUMSAMPLS, DAT ARRAY) 

« 

C  WALSH  TRANSFORM  SCALING 

* 

INTEGER  I 

•  INTEGER*2  NUMSAMPLS 

INTE6ER*4  DATARRAY(l) 

REAL  SCALEFACTOR 

•  * 

SCALEFACTOR  =  1.0  /NUMSAMPLS 
DO  10  I  ■  1,  NUMSAMPLS 

DATARRAY(l)  =  DATARRAYCI)  •  SCALEFACTOR 
10  CONTINUE 
RETURN 
END 


J 
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